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Executive  Summary 


Main  Result  and  Impact 

The  overall  aim  of  this  proposal  (FA9550-1 3-1 -0073)  was  to  address  the  cellular  mechanisms  of 
transcranlal  Direct  Current  Stimulation  (DCS).  Indeed,  the  scale  and  pace  on  tDCS  trials  In  humans 
(hundreds  published  per  year,  and  hundreds  ongoing)  has  vastly  outpaced  our  understanding  of  cellular 
mechanisms.  There  are  so  many  human  trials  because  various  Investigators  have  proposed  that  tDCS 
may  be  used  to  enhanced  wide  range  of  cognitive  functions  and  also  treat  a  wide  range  of  neuro- 
psychlatrlc  disease.  The  breadth  of  these  applications  Is  a/mosf  unbelievable  -  It  Is  based  on  the  general 
concept  that  tDCS  Is  paired  with  a  task  (some  kind  of  training)  and  that  tDCS  produces  an  enhancement 
In  only  that  task.  The  ability  of  tDCS  to  selectively  promote  performance  on  trained  tasks  Is  essential  but 
unaddressed  at  the  cellular  level.  The  combination  of  tDCS  with  training  implies  interactions  of  DCS  with 
ongoing  task-specific  activity  leading  to  task-specific  plasticity.  This  proposal  showed  that  DCS  can 
Induce  lasting  (plastic)  changes  selectivity  In  co-actIvated  synaptic  processes  and  provided  a  cellular 
substrate  for  these  changes.  This  proposal  therefore  provided  a  detailed  cellular-level  explanation  for  how 
tDCS  can  enhance  (theoretically)  any  task  It  Is  paired  with.  This  proposal  thus  not  only  provides  a  cellular 
explanation  for  tDCS,  It  also  showed  why  tDCS  could  be  applied  to  so  many  applications.  This  Is  a 
fundamental  advance  In  rationalizing  the  use  of  tDCS,  and  together  Is  a  new  theory  we  call  “functional 
targeting”. 


Report  Sections 

This  summary  report  Is  organized  on  chapters  that  approximated  papers  either  published  or  In 
preparation  that  resulted  from  this  proposal.  Related  to  the  overall  cellular  explant  Ion  for  tDCS 
summarized  above,  the  salient  results  of  each  chapter  are  simply  [Mapping  to  Alms  In  the  original 
proposal  Indicated]: 
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Section  1:  Identifies  a  new  cellular  target  of  tDCS.  Whereas  prior  work  has  focused  on  the  soma  of 
neurons,  we  show  the  axon  is  in  fact  more  sensitivity.  [Aim  1]. 

Section  2:  tDCS  is  diffuse  (current  going  to  much  of  the  brain)  and  weak  (low-intensity).  Whereas  prior 
work  tried  to  explain  away  this  “deficit”  we  in  fact  show  these  feature  of  tDCS  can  be  “benefits”  -  including 
making  the  effects  of  tDCS  specific  to  a  matched  training  task.  [Aim  1 ,  Aim  2]. 

Section  3:  tDCS  is  thought  to  boost  the  learning  of  tasks  or  therapy  applied  at  the  same  time.  We  provide 
a  cellular  mechanism  for  this.  Moreover,  we  show  that  thus  “boosting”  is  specific  to  the  trained  task.  [Aim 
2] 

Section  4:  tDCS  is  though  to  boost  learning  by  promoting  synaptic  plasticity.  We  show  the  mechanism  of 
this  plasticity.  [Aim  2] 
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Rahman  A,  Bikson  M.  et  al.  Electrical  stimulation  accelerates  and  boosts  the  capacity  for  synaptic 
learning.  In  preparation. 
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Bikson  M,  Truong  DQ,  Mourdoukoutas  AP,  Aboseria  M,  Khadka  N,  Adair  D,  Rahman  A.  Modeling 
sequence  and  quasi-uniform  assumption  in  computational  neurostimulation.  Prog  Brain  Res.  2015;222:1- 
23  doi:10.1016/bs.pbr.2015. 08.00 

Rahman  A,  Bikson  M.  Origins  of  specificity  during  tDCS:  anatomical,  activity-selective,  and  input-bias 
mechanisms  Frontiers  of  Human  Neuroscience  201 3;  doi  1 0.3389/fnhum.201 3.00688 

Reato  D,  Rahman  A,  Bikson  M,  Parra  LC.  Effects  of  weak  transcranial  Alternating  Current  Stimulation  on 
brain  activity  -  a  review  of  known  mechanisms  from  animal  studies.  Frontiers  of  Human  Neuroscience 
2013;  doi  10.3389/fnhum.201 3.00687 


Section  1:  Somatic  and  synaptic  terminal  effects  of  acute  direct  current 

stimuiation 

Abstract 

Transcranial  Direct  Current  Stimulation  (tDCS)  is  a  non-invasive  brain  stimulation  technique  to  modulate 
cortical  excitability.  Although  increased/decreased  excitability  under  the  anode/cathode  electrode  is 
nominally  associated  with  membrane  depolarization/hyperpolarization,  which  cellular  compartments 
(somas,  dendrites,  axons  and  their  terminals)  mediate  changes  in  cortical  excitability  remains 
unaddressed.  Here  we  consider  the  acute  effects  of  direct  current  stimulation  on  excitatory  synaptic 
efficacy.  Using  multi-scale  computational  models  and  rat  cortical  brain  slices  we  show:  1)  Typical  tDCS 
montages  produce  predominantly  tangential  (relative  to  the  cortical  surface)  direction  currents  (4-12  times 
radial  direction  currents),  even  directly  under  electrodes.  2)  Radial  current  flow  (parallel  to  the 
somatodendritic  axis)  modulates  synaptic  efficacy  consistent  with  somatic  polarization,  with  depolarization 
facilitating  synaptic  efficacy.  3)  Tangential  current  flow  (perpendicular  to  the  somatodendritic  axis) 
modulates  synaptic  efficacy  acutely  (during  stimulation)  in  an  afferent  pathway  specific  manner  that  is 
consistent  with  terminal  polarization,  with  hyperpolarization  facilitating  synaptic  efficacy.  4)  Maximal 
polarization  during  uniform  direct  current  stimulation  is  expected  at  distal  (branch  length  is  >3  times  the 
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membrane  length  constant)  synaptic  terminals,  independent  of  and  2-3  times  more  susceptible  than 
pyramidal  neuron  somas.  We  conclude  that  during  acute  direct  current  stimulation  the  cellular  targets 
responsible  for  modulation  of  synaptic  efficacy  are  concurrently  somata  and  axon  terminals,  with  the 
direction  of  cortical  current  flow  determining  the  relative  influence. 


Introduction 

During  tDCS,  current  flow  (~1  mA)  from  an  anode  to  a  cathode  electrode  generates  weak  electric  fields 
(<1  V/m)  across  the  cortex  (Datta  et  al.,  2009;  Salvador  et  al.,  2010).  tDCS  modulates  cortical  excitability 
in  the  primary  motor  cortex  (Nitsche  and  Paulus,  2000,  2001;  Antal  et  al.,  2004),  with  anodal  stimulation 
enhancing  and  cathodal  stimulation  diminishing  corticospinal  excitability,  as  measured  by  motor  evoked 
potentials  (MEPs)  elicited  by  transcranial  magnetic  stimulation  (TMS)  (Nitsche  et  al.,  2005).  Similarly,  in 
animal  models  of  DCS  spontaneous  and  evoked  cortical  potentials  were  acutely  facilitated  under  the 
anode  and  inhibited  under  the  cathode  (Creutzfeldt  et  al.,  1962;  Bindman  et  al.,  1964;  Purpura  and 
McMurtry,  1965b).  The  acute  changes  in  synaptic  efficacy  by  DCS  may  translate  to  enduring  effects 
(short  or  long-term  plasticity)  lasting  over  one  hour  after  stimulation,  dependent  on  the  duration  of 
stimulation  (typically  minutes)  (Bindman  et  al.,  1962,  1964;  Gartside,  1968)  and  the  nature  of  ongoing 
(synaptic)  activity  (Fritsch  et  al.,  2010;  Ranieri  et  al.,  2012).  Additionally,  the  acute  effects  of  DCS  are  not 
homogeneous  since  the  cellular  effects  of  stimulation  depend  on  neuronal  morphology,  stimulation 
intensity,  neuronal  orientation  relative  to  the  induced  electric  field,  and  on  the  nature  of  the 
spontaneous/evoked  activity(Tranchina  and  Nicholson,  1986;  Chan  et  al.,  1988;  Andreasen  and 
Nedergaard,  1996;  Bikson  et  al.,  2004;  Joucia  and  Yvert,  2009;  Radman  et  al.,  2009b).  Here,  we  consider 
if  characterizing  the  cellular  targets  of  DCS  may  help  reconcile  acute  neuromodulation  patterns  in  a  single 
framework.  We  specifically  focus  on  the  role  of  acute  cortical  DCS  on  presynaptic  (afferent  axon)  versus 
postsynaptic  (soma/dendrites)  cellular  compartments  in  modulating  synaptic  efficacy  (Jefferys,  1981; 
Bikson  et  al.,  2004;  Fritsch  et  al.,  2010;  Kabakov  et  al.,  2012;  Ranieri  et  al.,  2012). 
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Neuronal  excitability  in  resting  neurons  is  modulated  by  subthreshold  DCS  through  cell  membrane 
polarization  (<  1  mV  polarization  per  V/m)  (Radman  etal.,  2009).  While,  the  increase  in  excitability  under 
the  anode  is  attributed  to  “membrane  depolarization”  and  the  decrease  in  excitability  under  the  cathode  is 
attributed  to  “membrane  hyperpolarization”  (Bindman  et  al.,  1962,  1964;  Purpura  and  McMurtry,  1965b), 
in  fact,  during  DCS  there  are  an  equal  number  of  cellular  elements  that  are  hyperpolarized  or  depolarized 
in  any  given  brain  region  (Joucia  &  Yvert,  2009),  including  directly  under  the  anode  and  cathode.  For 
example,  during  DCS  of  cortical  pyramidal  neurons,  somatic  depolarization  is  associated  with  concurrent 
apical  dendritic  hyperpolarization,  and  somatic  hyperpolarization  is  associated  with  apical  dendritic 
depolarization  (Chan  and  Nicholson,  1986;  Tranchina  and  Nicholson,  1986;  Chan  et  al.,  1988;  Andreasen 
and  Nedergaard,  1996;  Radman  et  al.,  2009b).  Additionally,  afferent  axons  and  their  terminals  are  also 
polarized  during  DCS  with  varied  polarity  depending  on  morphology  (BeMent  and  Ranck,  1969;  Ranck, 
1975;  Tranchina  and  Nicholson,  1986;  Arlotti  et  al.,  2012;  Kabakov  et  al.,  2012;  Marquez-Ruiz  et  al., 
2012).  Finally,  membrane  compartment-specific  polarization  affects  the  capacity  of  presynaptic  input  to 
influence  postsynaptic  output  (synaptic  efficacy)  in  both  acute  and  long-lasting  effects  of  DCS  (Bikson  et 
al.,  2004;  Fritsch  et  al.,  2010;  Kabakov  et  al.,  2012;  Marquez-Ruiz  et  al.,  2012).  The  diversity  of  cellular 
targets  of  DCS  raises  a  number  of  mechanistic  questions.  Which  membrane  compartments  (somas, 
dendrites,  or  axons/terminals)  polarized  during  DCS  contributes  to  changes  in  excitability,  including 
synaptic  efficacy?  Does  depolarization/hyperpolarization  of  the  implicated  compartments  in  fact  correlate 
with  facilitation/inhibition  of  synaptic  efficacy? 

Classical  animal  research  has  implicated  somatic  polarization  during  transcortical  DCS  (Terzuolo  and 
Bullock,  1956;  Purpura  and  McMurtry,  1965a,  b).  Accordingly,  local  field  potentials  and  intracellular 
recordings  suggest  that  facilitation/inhibition  of  spontaneous  activity  (Frohlich  &  McCormick,  2010;  Reato 
et  al.,  2010)  and  synaptic  efficacy  during  DCS  may  be  associated  with  depolarization/hyperpolarization  of 
the  somatic  membrane  potential  (Jefferys,  1981;  Bikson  et  al.,  2004;  Radman  et  al.,  2009).  Flowever, 
afferent  axons,  including  cortico cortical  and  thalamocortical  connections,  which  are  involved  in  synaptic 
processing,  might  additionally  contribute  to  the  effects  of  DCS  (Purpura  and  McMurtry,  1965a,  b;  Bikson 
et  al.,  2004;  Kabakov  et  al.,  2012;  Marquez-Ruiz  et  al.,  2012).  For  example,  modulation  of  presynaptic 
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transmitter  release  at  thalamocortical  sensory  afferents  has  been  attributed  to  the  effects  of  DCS  in  vivo 
(Marquez-Ruiz  et  al.,  2012).  Whether  hyperpolarization  or  depolarization  of  axon  terminals  enhances 
efficacy  may  vary  across  synaptic  pathways.  Finally,  DCS  of  the  apical  dendrites  of  pyramidal  neurons, 
which  polarize  opposite  to  somas,  may  further  influence  synaptic  processing  (Bikson  et  al.,  2004).  Thus,  it 
remains  unclear,  which  compartments  influence  modulation  by  DCS  (somas,  dendrites,  or 
axons/terminals)  and  whether  depolarization  or  hyperpolarization  is  associated  with  enhanced  synaptic 
efficacy.  Understanding  the  cellular  targets  of  tDCS  is  considered  pre-requisite  to  rational  electrotherapy 
design  and  optimization  (Bikson  et  al.,  2012b). 

The  cellular  targets  of  TES  remain  unclear.  Afferent  axons,  somas,  and  dendrites  are  all  susceptible  to 
polarization  depending  on  the  directionality  of  the  induced  electric  field  and  the  orientation  of  the  neuronal 
compartment.  The  well-established  approach  to  characterize  the  cellular  effects  of  electrical  stimulation  is 
to  determine  the  electric  field  produced  in  the  brain  target  and  then  reproduce  the  electric  field  in  animal 
models  or  neuron  compartment  models;  for  transcranial  stimulation,  including  with  DCS,  the  quasi¬ 
uniform  assumption  is  widely  applied  (Bikson  et  al.,  2012b).  We  extended  previous  MRI-derived  models 
of  tDCS  generated  brain  current  flow  (Datta  et  al.,  2009;  Salvador  et  al.,  2010)  to  quantify  the 
directionality  of  the  induced  electric  field  in  the  cortex  -  specifically  the  relative  radial  and  tangential 
components  of  the  cortical  electric  field.  In  the  second  part  of  this  study,  we  aimed  to  characterize  the  role 
of  uniform  electric  fields  on  modulating  the  excitability  of  extracellular  field  potentials  using  brain  slices, 
which  facilitate  control  of  electric-field  direction  (radial  and  tangential  to  the  somatodendritic  axis)  and 
isolated  activation  of  pathway-specific  synaptic  activity.  Finally,  compartment  models  of  morphologically 
reconstructed  cortical  pyramidal  neurons  were  used  to  predict  which  neuronal  elements  are  polarized  by 
uniform  DC  fields. 

We  report  that  during  tDCS  both  radial  and  tangential  (relative  to  the  cortical  surface)  direction  currents 
are  induced,  however,  tangential  currents  are  prevalent  across  the  cortex.  Fields  radial  to  the  cortical 
surface  modulated  synaptic  efficacy  independent  of  synaptic  pathway  and  consistent  with  somatic 
polarization,  with  depolarization/hyperpolarization  facilitating/inhibiting  synaptic  efficacy,  while  tangential 
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fields  modulated  synaptic  efficacy  in  a  pathway-specific  manner,  such  that  axon  terminal 
hyperpolarization/depolarization  was  associated  with  facilitation/inhibition  of  synaptic  efficacy.  The  acute 
effects  of  DCS  on  synaptic  efficacy  may  be  reconciled  through  (concurrent)  contributions  from  these  two 
cellular  targets. 

The  overall  goal  of  this  study  was  to  systematically  characterize  the  cellular  target  of  action  for  DC 
stimulation.  Our  approach  was  to  first  determine  the  prevalence  of  radial  and  tangential  EFs  in  the  cortex 
using  MRI-derived  computational  models.  Secondly,  effects  of  radial  and  tangential  DC  fields  on  synaptic 
efficacy  were  characterized  in  cortical  brain  slices  in  four  distinct  synaptic  pathways  corresponding  to 
different  orientations  of  afferent  axons.  Synaptic  efficacy  was  shown  to  be  pathway-specific  and  synaptic 
terminal  polarization  by  tangential  fields  can  lead  to  facilitation/inhibition  independent  of  somatic 
polarization.  These  findings  motivated  us  to  determine  a  general  rule  to  describe  sensitivity  of  axon 
terminals  to  polarization  in  the  form  of  a  coupling  constant. 


Methods 

MRI-derived  electric  field  models 

Complete  methods  for  MRI-derived  head  models  of  electric  field  (EF)  distributions  are  detailed  elsewhere 
(Datta  et  al.,  2009).  Briefly,  EF  distributions  in  the  brain  were  computed  using  individualized  head  models 
created  from  the  T1 -weighted  MRI  scans  of  a  healthy  individual.  The  head  was  segmented  into 
compartments  representing  gray  matter,  white  matter,  cerebrospinal  fluid  (CSF),  skull,  scalp,  eye  region, 
muscle,  air,  and  blood  vessels  using  a  combination  of  tools  from  the  FSL  and  Simpleware  (Exeter,  UK). 
The  finite  element  (FE)  mesh  generated  from  the  segmentation  masks  was  exported  to  COMSOL 
Multiphysics  4.2  (Burlington,  MA)  for  computation  of  electric  fields. 

Electric  field  magnitudes  were  calculated  for  the  conventional  tDCS  montage  over  M1-SO  (one  anode 
over  the  primary  motor  cortex  and  one  cathode  over  the  supraorbital  region,  5x5  cm  sponge  electrodes. 
Fig.  1A1,  2A)  and  the  Fligh-Definition  tDCS  (FID-tDCS)  montage  (center  anode  over  the  primary  motor 
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cortex  and  four  surround  cathode  electrodes,  1  cm  electrode  diam.  with  6  cm  separation  between  center 
and  surround  electrodes)  (Datta  et  al.,  2009).  The  relative  magnitude  of  the  two  components  of  the  EF  ( 
=  radial  and  =  tangential)  are  considered  and  quantified  on  multiple  scales  (Fig.  10),  Including 

global  field  distributions  In  the  brain,  regionally  under/between  electrodes,  and  In  sub-regions  on  gyral 
crowns/walls.  We  use  the  ratio  of  tangential  to  radial  {E^IE^)  field  magnitudes  to  describe  the  relative 

magnitudes  In  each  region,  such  that  E^l E^>1  corresponds  to  greater  tangential  fields  on  average  and 


E^/E^<1  corresponds  to  greater  radial  fields  on  average  (Fig.  9A2,  A3).  The  E^IE^  metric  Is 


represented  In  Fig.  1A1  with  a  schematic  representation  of  the  voltage  distribution  overlaid  on  each 


region  of  Interest  along  a  cortical  gyrus. 
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Figure  1.  Multi-scale  methods  and  outcome  measures  of  uniform  electric  field  directionality  and 
effects.  A1,  GyrI-precIse  finite  element  models  of  current  flow  during  tDCS  Indicate  a  uniform 
voltage  gradient  In  cortical  grey  matter  (GM)  directly  under  the  anode.  A2,  The  Induced  electric 
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field  in  the  cortex  can  be  decomposed  into  a  radial  component  {E^)  that  is  parallel  to  the 
somatodendritic  axis  and  a  tangential  component  {E^)  that  is  orthogonal  to  the  somatodendritic 

axis.  A3,  We  quantified  the  relative  occurrence  of  radial  and  tangential  fields  in  cortical  GM 
expressed  as  the  ratio  of  the  average  of  the  field  magnitude  in  the  tangential  direction  to  the 
average  of  the  field  magnitude  in  the  radial  direction  {E^/  E^).  B1-2,  The  brain  slice  preparation 

was  used  to  study  the  change  in  synaptic  efficacy  during  a  uniform  radial  or  tangential  field  by 
recording  evoked  field  potentials.  The  voltage  gradient  between  parallel  Ag/AgCI  wires  is 
superimposed  on  a  schematic  of  a  sagittal  slice  of  the  rat  primary  motor  cortex.  From  the 
macroscopic  to  the  mesoscopic  scale  we  can  approximate  a  uniform  electric  field  along  the  length 
of  a  neuron  (compare  voltage  gradients  in  A1  and  B1).  B3,  The  field  EPSP  provides  a  measure  of 
synaptic  efficacy  through  facilitation  or  inhibition  of  the  response  amplitude.  C1,  Compartment 
model  simulations  of  morphologically  reconstructed  neocortical  pyramidal  neurons  were  used  to 
provide  a  description  of  axon  terminal  polarization  in  a  uniform  electric  field.  C2,  The  polarization 
profile  of  a  layer  5  pyramidal  neuron  in  a  radially  directed  uniform  electric  field  indicates  soma 
depolarization  (red)  corresponds  to  apical  dendrite  hyperpolarization  (blue).  Layer  2/3  neurons 
have  a  more  complex  polarization  profile  with  long  processes  reaching  maximal  depolarization 
independent  of  the  neuronal  body.  C3,  Neurons  in  a  uniform  electric  field  directed  tangential  to 
the  somatodendritic  axis  preferentially  affects  processes  that  are  oriented  along  the  tangential 
field. 


Electrophysiology 

Brain  slices  including  a  part  of  the  motor  cortex  were  prepared  from  male  young  adult  Wistar  rats  aged  3 
to  6  weeks  old,  which  were  deeply  anesthetized  with  ketamine  (7.4  mg/kg)  and  xylazine  (0.7  mg/kg) 
applied  intraperitoneally  (i.p.)  and  sacrificed  by  cervical  dislocation.  The  brain  was  removed  and 
immersed  in  chilled  (2-6  °C)  artificial  cerebrospinal  fluid  (ACSF)  containing  (in  mM):  126  NaCI,  3  KCI,  1.25 
NaFl2P04,  2  MgS04,  2  CaCI,  24  NaFICOs  and  10  D-glucose,  bubbled  with  a  mixture  of  95%  02-5%  CO2). 
Parasagittal  slices  (450  pm  thick)  were  cut  at  the  distance  of  2. 0-3.0  mm  from  the  brain  midline  using  a 
vibrating  microtome  and  transferred  to  a  holding  chamber  for  at  least  1  h  in  ambient  temperature.  Slices 
were  then  transferred  to  a  fluid-gas  interface  chamber  perfused  with  warmed  ACSF  (30.0±0.5  °C)  at  1 .9 
mL/min.  The  humidified  atmosphere  over  the  slices  was  saturated  with  a  mixture  of  95%  02-5%  C02. 
Recordings  started  2-3  h  after  dissection. 


To  investigate  the  role  of  distinct  cortical  pathways  on  modulation  of  synaptic  efficacy  the  recording 
electrode  (a  glass  micropipette  filled  with  0.25  M  NaCI,  resistance  1-8  MQ)  was  placed  in  either  layer  ll/lll 
or  layer  V  of  the  rat  primary  motor  cortex  (Fig.  11  A).  Neocortical  horizontal  pathways  (within  a  cortical 
layer)  and  vertical  pathways  (across  lamina)  represent  functionally  distinct  connections  that  are  well 
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characterized  in  the  rat  neocortex  (Aroniadou  and  Keiier,  1993;  Keiier,  1993).  For  instance,  horizontai 
connections  integrate  information  across  sensory  maps  in  the  motor  cortex  by  connecting  reiated 
functionai  coiumns;  however,  extrinsic  connections  from  distant  corticai  and  subcorticai  regions  have 
dominantiy  verticai  projections  (Hess  and  Donoghue,  1994;  Sasaki  et  ai.,  2011).  The  synaptic  circuits 
representing  intrinsic  connections  between  popuiations  of  neurons  within  the  motor  cortex  was  previousiy 
identified  in  the  iiterature  using  current  source  density  distributions,  morphoiogicai  anaiysis,  and  through 
studies  of  functionai  connectivity,  inciuding  iong-term  potentiation  and  depression  experiments 
(Aroniadou  and  Keiier,  1993;  Keiier,  1993;  Aroniadou  and  Keiier,  1995;  Castro-Aiamancos  et  ai.,  1995; 
Hess  et  ai.,  1996;  Riouit-Pedotti  et  ai.,  1998).  Fig.  3B  summarizes  the  identified  pathways  reievant  for  this 
study,  in  addition  to  probing  the  pathway-specific  effects  of  eiectric  fieids  on  corticai  fieid  potentiais,  we 
use  the  synaptic  organization  of  the  cortex  to  test  the  effectiveness  of  radiai  and  tangentiai  fieids  in 
moduiating  synaptic  efficacy  in  horizontai  and  verticai  pathways. 

Orthodromic  stimuiation  targeting  four  distinct  synaptic  pathways  previousiy  identified  in  the  rat  primary 
motor  cortex  was  appiied  with  a  bipoiar  piatinum/stainiess  steei  stimuiating  eiectrode  piaced  either  400- 
800  pm  (stimuiating  eiectrodes  SI  and  S2)  or  1100-1300  pm  (stimuiating  eiectrodes  S3  and  S4)  beiow 
the  piai  surface  to  activate  fibres  running  within  iayer  ii/iii  or  iayer  V,  respectiveiy  (Fig.  11  A).  The 
stimuiating  eiectrodes  were  piaced  either  iateraiiy  anterior  or  posterior  to  the  recording  eiectrode  targeting 
horizontai  afferent  synaptic  connections  in  either  of  the  superficiai  (SI  and  S2)  or  deep  iayers  (S3)  (Fig. 
11  A).  Verticai  connections  were  stimuiated  with  a  bipoiar  eiectrode  (S4)  in  deep  iayer  V  and  a  recording 
eiectrode  in  iayer  ii/iii.  Fieid  excitatory  postsynaptic  potentiais  (fEPSPs),  which  provides  a  measure  of 
iocaiized  extraceiiuiar  currents  generated  from  a  popuiation  of  pyramidai  neurons  in  response  to  synaptic 
or  antidromic  activation,  were  evoked  with  constant-current  puises  (0.2  ms)  deiivered  once  per  minute. 
The  test  stimuius  intensity  (30-200  pA)  was  adjusted  to  evoke  haif-maximai  responses  based  on  input- 
output  curves  (Fig.  12A);  no  consistent  reiationship  was  found  between  fEPSP  deiay  and  stimuiation 
ampiitude  or  fEPSP  peak.  This  reiativeiy  weak  stimuiation  did  not  evoke  popuiation  spikes,  and  usuaiiy 
eiicited  a  fEPSP  having  a  singie  peak  in  84  of  92  siices  tested.  Responses  were  ampiified,  iow  pass 
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filtered  (1000  Hz  cutoff,  Warner  Instruments,  CT,  USA),  acquired  at  a  10-kHz  sampling  rate  (1401 
interface,  CED,  UK)  and  analyzed  on-  and  offline  (Signal  3  software,  CED,  UK). 

The  stimulation  protocol  was  designed  to  measure  the  acute  (during  field)  changes  in  fEPSPs  evoked 
during  an  applied  electric  field.  fEPSP  amplitudes  were  used  as  a  measure  of  the  change  in  synaptic 
efficacy  during  the  field,  such  that  an  increase  in  the  amplitude  indicates  facilitation  and  a  decrease 
indicates  inhibition.  Uniform  fields  (±8  V/m)  were  generated  by  passing  current  (D/A  driven  analog  current 
follower,  A-M  Systems,  WA,  USA)  between  two  large  parallel  Ag-AgCI  wires  positioned  in  the  bath  across 
the  slice  for  1  sec  starting  0.5  sec  before  bipolar  stimulation.  The  field  intensity  was  chosen  based  on  pilot 
experiments  to  produce  approximately  a  5%  change  in  fEPSP  amplitude  such  that  a  statistically 
significant  change  could  be  observed  within  slices  with  a  practical  sample  size  (Fig.  12B). 

The  quasi-uniform  EF  assumption  (Bikson  et  al.,  2012b),  which  considers  the  applied  EF  amplitude  may 
be  uniform  on  the  scale  of  the  membrane  length  constant,  allowed  us  to  evaluate  synaptic  efficacy  in  a 
population  of  neurons  in  a  uniform  EF  in  vitro  as  well  as  in  single  neuron  models  (compare  Fig.  1A1,  B1 
and  Cl).  Specifically,  tDCS  produces  approximate  uniform  fields  across  select  regions  of  grey  matter 
(Fig.  12A,  linear  voltage  gradients)  that  can  be  reproduced  in  vitro  (Fig.  12B)  and  simulated  in  silico  (Fig. 
12C).  Consideration  of  uniform  electric  fields  allows  analysis  in  this  study  across  the  macroscopic  (brain 
current  flow  using  FEM),  mesoscopic  (synaptic  efficacy  in  brain  slice),  and  microscopic  scales  (neuron 
compartment  polarization).  Indeed,  representation  of  electric  fields  is  ubiquitous  in  FEM  computational 
studies  (Miranda  et  al.,  2006;  Datta  et  al.,  2009;  Salvador  et  al.,  2010;  Datta  et  al.,  2011),  and  uniform 
electric  fields  are  used  in  mechanistic  studies  of  tDCS  (Radman  et  al.,  2009a;  Anastassiou  et  al.,  2010; 
Fritsch  et  al.,  2010;  Arlotti  et  al.,  2012;  Kabakov  et  al.,  2012;  Ranieri  et  al.,  2012).  To  our  knowledge  this 
is  the  first  study  to  integrate  analysis  across  these  three  scales. 

All  the  statistical  data  are  expressed  as  mean  ±  SD,  unless  stated  otherwise.  The  statistical  difference 
between  groups  (critical  value=0.05)  was  estimated  using  the  Student’s  t-test  considering  normally 
distributed  fEPSP  amplitudes  (Lilliefors  test  for  normality). 
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Compartment  model  simulation 

Numerical  and  analytical  solutions  of  3D  reconstructed  neurons  and  simplified  cable  models,  respectively, 
were  used  to  infer  a  general  rule  of  steady-state  axon  terminal  polarization  under  a  uniform  electric  field 
based  on  axonal  morphology.  With  thousands  of  axonal  afferents,  with  diverse  morphologies  for  any 
given  cortical  neuron,  a  range  of  axon  terminal  polarization  values  are  expected;  but  our  goal  was  to 
approximate  the  maximum  polarization  expected  in  the  most  sensitive  axon  terminals.  Specifically,  we 
aimed  to  characterize  the  sensitivity  of  axon  terminals  by  their  coupling  constant  (which  linearly  relates 
steady-state  membrane  polarization  to  electric  field);  this  then  allows  comparison  of  maximal  axon 
terminal  sensitivity  with  previously  characterized  coupling  constants  for  neuron  somas  (Radman  et  al., 
2009a). 

Numerical  solutions  for  realistic  neurons.  Complex  neuronal  geometries  require  numerically  solving  the 
cable  equation  (Basser  and  Roth,  2000).  Eight  morphologically  reconstructed  layer  ll/lll  and  layer  V  rat 
somatosensory  pyramidal  neurons  from  the  NeuroMorpho.org  database  (ID:  NMO_0116,  NMO_00348, 
NMO_00410,  NMO_00417,  NMO_01134,  NMO_01132,  NMO_01135,  NMO_01112)  were  used  to 
approximate  membrane  polarization  in  a  uniform  EF  using  Neuron  (version  7.2)  (Wang  et  al.,  2002; 
Ascoli,  2006).  The  built-in  function  extracellular(x)  was  used  to  apply  1  V/m  uniform  electric  fields  in  the 
radial  direction  (parallel  to  the  somatodendritic  axis.  Fig.  1C1,  C2),  tangential  direction  (perpendicular  to 
the  somatodendritic  axis.  Fig.  1C3),  and  in  optimal  polarization  angles  along  individual  axons  to  produce 
the  maximum  terminal  polarization  in  that  axon.  The  membrane  resistance  was  70  kOcm"^  and  the 
intracellular  axial  resistance  R.  was  155  Qcm  (Markram  et  al.,  1998).  The  median  axonal  length  constant 

{2.=  I - ^  ,  where  r  is  the  radius  p  is  the  membrane  resistivity  and  p.  is  the  intracellular 

V2xp. 

resistivity)  was  0.56  mm  (0.05  mm  interquartile  range).  In  all  compartment  neuron  models,  passive 
resistive  parameters  were  used  to  approximate  steady  state  membrane  polarization. 
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Branched  semi-infinite  axon.  Next,  we  considered  a  bent  semi-infinite  axon  in  a  uniform  EF  to  establish  a 


relationship  between  terminal  polarization  and  fibre  morphology,  extending  previous  derivations 
(Tranchina  and  Nicholson,  1986;  Basser  and  Roth,  2000;  Miranda  et  al.,  2007b).  We  previously  derived 
the  analytical  steady-state  solution  to  the  classic  cable  equation  for  a  passive  terminating  straight  branch 
of  an  axon  (Eq.  1,  for  a  complete  derivation  see  (Arlotti  et  al.,  2012). 

Equation  1.  V)  =  £’/lcos(0)tanh(i'//l)-l-yo - ^ - 

cosh(i/A) 

The  terminal  polarization  varies  with  the  physical  length  £  of  the  final  segment  and  depends  on  the  angle 
6  of  the  electric  field  E  relative  to  the  final  branch  and  the  membrane  polarization  at  the  bend  Vq  .  The 

relative  contribution  of  EX  vs.  depends  on  the  electrotonic  (dimensionless)  length  IjX  . 

Straight  semi-infinite  axon.  A  simple  analytical  solution  of  polarization  of  a  straight  semi-infinite  axon  in  a 
uniform  EF  is  well  known  (and  is  in  fact  a  special  case  of  the  axon  branch  with  infinite  length);  Equation  2 
relates  polarization  at  the  fibre  terminal  with  the  membrane  space  constant  X  (Flause,  1975).  The 
terminal  membrane  polarization  is  a  function  of  the  uniform  extracellular  EF  relative  to  the  axon,  such  that 
6  =  0  corresponds  to  an  EF  along  the  longitudinal  axis  of  the  fibre.  Interestingly,  in  this  formulation  X  is 
analogous  to  the  coupling  constant,  and  we  intend  to  establish  if  and  when  this  approximation  is  valid  for 
realistic  axon  morphologies. 


Equation  2. 


^  »  A  ^  y,  =  EXcosO 


Results 
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Quantitative  analysis  of  electric  field  directionality  in  finite  element  models 

Consistent  with  previous  computational  studies,  tDCS  and  HD-tDCS  are  predicted  to  generate  weak  EFs 
(<1  V/m  per  mA  applied  current)  in  the  cortex,  with  a  complex  spatial  pattern  determined  by  the  electrode 
montage  and  the  details  of  brain  anatomy,  notably  cortical  surface  idiosyncrasies  (Datta  et  al.,  2009). 
Under  the  assumption  that  the  voltage  gradient  in  each  region  of  the  grey  matter  was  locally  linear,  we 
summarized  the  (uniform)  EF  in  each  region  (Fig  12A,  Fig  13).  Cortical  EFs  generated  during  stimulation 
have  components  both  radial  {E^  ,  normal)  and  tangential  {E^, ,  parallel)  to  the  cortical  surface  (Fig.  9B). 

Radial  fields  that  are  oriented  along  the  somatodendritic  axis  of  cortical  pyramidal  neurons  produce 
somatic  polarization.  Tangential  fields  are  transverse  to  the  somatodendritic  axis  of  cortical  pyramids,  and 
will  polarize  horizontally  directed  corticocortical  axons.  Fligh-resolution  modeling  shows  that  current 
flowing  across  grey  matter  will  have  both  radial  and  tangential  components  in  conventional  (Fig.  10A,  top) 
and  Fligh-Definition  (Fig.  10A,  bottom)  tDCS.  We  previously  reported  interindividual  variations  of  peak 
cortical  electric  field  of  -two  fold  across  anatomically  normal  adult  heads  (Datta  et  al.,  2012)  but  focus 
here  on  relative  directional  distribution.  The  global  prevalence  of  tangential  fields  motivated  further 
analysis  of  field  directionality  on  additional  levels  including  regions  under/between  electrodes  and  sub- 
regions  of  gyral  crown/walls. 

The  distribution  of  radial  and  tangential  components  (Fig.  10B;  density  plots  show  relative  occurrence) 
under  the  anode  {EJE^  tDCS:  7.9,  HD-tDCS:  4.8),  cathode  {EJE^  tDCS:  7.0,  HD-tDCS:  8.0),  and 

between  electrodes  {E^,IE^  tDCS:  11.7,  HD-tDCS:  7.6)  indicate  most  elements  have  both  field 

components  but  with  (significantly)  stronger  tangential  than  radial  field  (Fig.  10B,  inset  histograms). 
Interestingly,  the  few  isolated  highest  EFs  are  radial  (Fig.  10B,  axis  histograms). 

Cortical  folding  greatly  influenced  the  EF  distribution  under  the  electrodes  (Fig.  IOC).  For  example,  on 
the  scale  of  a  single  gyrus  in  the  primary  motor  strip,  the  gyral  crown  has  a  dominantly  tangential 
component  {E^IE^  tDCS:  16.6,  HD-tDCS:  7.3)  and  the  strongest  fields  are  tangentially  oriented.  In 

contrast,  motor  strip  gyral  walls  have  a  relatively  equal  distribution  of  field  components  using  the 
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conventional  tDCS  electrode  montage  and  slightly  stronger  radial  components  in  HD-tDCS  {E^/E^ 

tDCS:  1.1,  HD-tDCS:  0.33).  The  maximal  EF  in  gyral  walls  is  typically  radial  to  the  cortical  surface,  but  a 
few  elements  have  negligible  tangential  components  as  well. 

In  all  regions,  under  and  between  electrodes,  there  is  a  higher  density  of  elements  with  stronger 
tangential  than  radial  fields,  and  few  elements  had  purely  radially  oriented  fields.  Given  the  prevalence  of 
tangential  EFs,  we  next  considered  field  direction  sensitive  changes  in  synaptic  efficacy  in  cortical 
pathways  of  the  rat  primary  motor  cortex. 
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Figure  2.  Forward  model  of  tDCS  and  FID-tDCS  quantifying  electric  field  direction  metrics.  During 
tDCS,  current  may  be  dominantly  tangential  (along  the  cortical  surface)  rather  than  radial,  even  in 
brain  regions  directly  under  the  electrodes.  A,  MRI-derived  finite  element  models  of  current 
distribution  in  a  gyri-precise  head  model  are  used  to  quantify  the  relative  occurrence  of  radial 
(normal  to  the  cortical  surface)  and  tangential  (along  the  cortical  surface)  components  of  the 
electric  field.  Both  conventional  (top)  and  high  definition  (FID,  bottom)  tDCS  montages  produce 
radial  {E^ ,  normal  to  the  cortical  surface)  and  tangential  current  {E^ ,  along  the  cortical  surface) 

indicated  by  the  global  electric  field  distribution  (V/m)  across  the  head.  In  the  FID-tDCS  montage, 
current  is  focalized  within  the  ring  configuration  (inset)  with  radial  currents  under  the  center 
electrode  and  tangential  currents  between  the  surround  electrodes.  Qualitative  comparison  of  the 
electric  field  components  indicate  greater  radial  field  magnitudes  in  the  gyral  wall  and  greater 
tangential  field  magnitudes  in  the  gyral  crown  (compare  insets).  B,  Regionally,  the  distribution  of 
field  component  magnitudes  indicate  prevailing  tangential  currents  under  the  anode,  cathode  and 
between  electrodes  as  described  by  the  ratio  of  tangential  to  radial  field  magnitude  {E^/  E^  ratio. 
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see  methods).  However,  most  elements  have  both  radial  and  tangential  components,  and  the 
isolated  highest  electric  fields  are  radial.  The  tangential  and  radial  component  for  individual 
elements  is  shown  for  each  sub-region  in  false  color  density  plots,  which  show  relative 
occurrence  (relative  density  from  absent  (green)  to  maximal  (red)).  Axis  histograms  show  relative 
distribution  of  elements  with  a  given  tangential  or  radial  component  electric  field.  Inset  histograms 
describe  the  distribution  of  the  %  of  elements  in  a  region  as  a  function  of  the  normalized 
component  magnitude  (such  that  1  indicates  elements  with  dominant  radial  or  tangential 
component).  C,  Cortical  folding  further  influences  the  distribution  of  the  electric  field,  therefore, 
sub-regional  field  component  distributions  are  indicated  for  a  gyral  crown  and  wall.  Tangential 
fields  are  dominant  in  magnitude  in  the  gyral  crown  but  are  weaker  in  the  walls  where  radial 
magnitudes  are  stronger,  as  observed  in  A. 


Effects  of  parallel  and  perpendicular  electric  fields  on  field  potentials 

The  primary  motor  cortex  is  characterized  by  functionally  distinct  afferent  synaptic  pathways  with  specific 
axonal  morphologies  and  orientations  (Aroniadou  and  Keller,  1993;  Keller,  1993).  Using  the  rat  primary 
motor  cortex  slice  preparation  we  tested  the  acute  effects  of  uniform  EFs  on  synaptic  efficacy  in  four 
distinct  cortical  pathways  for  electric  fields  applied  parallel  to  the  somatodendritic  axis  (radial  fields)  and 
perpendicular  to  the  somatodendritic  axis  (tangential  fields).  Field  excitatory  postsynaptic  potentials  were 
monitored  in  layer  ll/lll  or  layer  V  in  response  to  activity  evoked  by  stimulation  of  posterior  or  anterior 
afferent  synaptic  pathways  in  layer  ll/lll  (stimulating  electrodes  SI  and  S2,  respectively)  and  horizontal  or 
vertical  afferents  in  layer  V  (stimulating  electrodes  S3  and  S4,  respectively)  (Fig.  11  A,  B).  In  addition  to 
monosynaptic  fEPSPs  (4-7  ms  peak),  a  fiber  volley  was  observed  in  some  slices  that  was  non- 
synaptically  mediated  as  indicated  by  the  time  course  (<2  ms  peak)  and  insensitivity  to  bath  application  of 
the  non-NMDA  receptor  antagonist  6,7-dinitroquinoxaline-2,3-dione  (DNQX,  20  pM)  (Fig.  1 1 D).  We  tested 
the  hypothesis  that  synaptic  efficacy  increases  with  somatic  depolarization  and/or  axon  terminal 
hyperpolarization;  thus  in  a  defined  polarity,  field  direction,  and  pathway-specific  manner  (Fig  11C). 
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Figure  3.  Electrophysiology  of  direction-specific  uniform  DC  electric  fields  in  synaptic  pathways  of 
the  rat  motor  cortex.  A,  Schematic  of  electrophysiology  setup  where  uniform  extracellular  electric 
fields  were  generated  in  all  experiments  by  passing  constant  current  across  parallel  Ag-AgCI 
wires  positioned  in  the  bath  across  the  slice.  Activity  was  monitored  in  layer  ll/lll  (REF)  or  layer  V 
with  a  glass  microelectrode.  An  additional  field  electrode  (REF)  was  positioned  in  an  iso-potential 
to  remove  the  uniform  field  artifact.  Activity  was  evoked  with  a  bipolar  stimulating  electrode  (S1- 
S4)  positioned  500  pm  from  the  recording  electrode  in  either  layer  ll/lll  or  layer  V  targeting  one  of 
four  distinct  synaptic  pathways  corresponding  to  different  orientations  of  afferent  axons:  posterior 
horizontal  layer  ll/lll  (S1),  anterior  horizontal  layer  ll/lll  (S2),  posterior  horizontal  layer  V  (S3),  and 
vertical  layer  V  to  ll/lll  (S4).  B,  Diagram  summarizing  the  primary  synaptic  circuits  in  this  study. 
Line  thickness  and  diameter  of  the  filled  circles,  which  represent  synapses,  are  correlated  with 
the  strength  of  the  synaptic  input.  C,  Schematic  of  the  expected  polarization  in  distinct  synaptic 
pathways  exposed  to  radial  and  tangential  fields.  Somas,  dendrites,  axons  and  axon  terminals 
are  depolarized  (red),  hyperpolarized  (blue),  or  not  affected  (black)  by  DC  fields.  D,  Characteristic 
fEPSP  and  field  spike  waveforms  from  the  layer  V  pathway.  The  fEPSPs,  but  not  earlier  field- 
spike,  were  suppressed  by  the  non-NMDA  receptor  antagonist  DNQX. 
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Radial  EFs  (parallel  to  the  somatodendritic  axis  of  cortical  pyramidal  neurons)  typically  leads  to 
‘excitation’  for  positive  fields  (anode  proximal  to  pia)  and  ‘inhibition’  for  negative  fields  (cathode  proximal 
to  pia)  (Fig.  9B3).  In  three  of  the  afferent  pathways  tested  fEPSP  amplitudes  were  significantly  facilitated 
(Fig.  12C;  S1:  posterior  horizontal  layer  ll/lll,  5.2±3%,  p=0.031,  n=7;  S3:  horizontal  layer  V,  7.3±8%, 
p=0.02,  n=7;  S4:  vertical  layer  ll/lll,  9.4±4%  p=0.023,  n=11)  by  radial  positive  EFs  (+8  V/m),  with  no 
significant  effect  in  one  pathway  (S2:  anterior  horizontal  layer  ll/lll,  p=0.251 ,  n=6).  Radial  negative  fields  (- 
8  V/m)  reduced  responses  in  the  same  three  pathways  (S1:  -8±1%,  p=0.022,  n=7;  S3:  -6.4±6.9%, 
p=0.024,  n=7;  S4:  -10.2±1.5%,  p=0.01,  n=11)  with  no  significant  effect  in  the  anterior  horizontal  layer  ll/lll 
pathway  (S2:  p=0.98,  n=6).  A  change  in  fEPSP  timing  was  not  resolved.  Importantly,  grouping  across  all 
pathways  (e.g.  not  controlling  pathway  selectivity),  8  V/m  radial  positive  fields  significantly  facilitated 
(7±1%,  p=0.004)  and  radial  negative  fields  inhibited  (-6.1  ±4%,  p=0.04)  fEPSP  amplitudes  (Fig.  12C). 

tDCS  also  induces  tangential  fields  that  are  not  expected  to  polarize  neurons  with  a  somatodendritic  axis 
radial  to  the  cortical  surface  (Fig.  12C)  (Bikson  et  al.,  2006).  Yet,  we  found  tangential  positive  fields 
significantly  modulated  three  of  four  pathways  tested;  moreover,  the  direction  of  change  in  excitability  was 
pathway-specific.  Tangential  positive  fields  (+8  V/m,  oriented  posterior  to  anterior)  reduced  fEPSP 
responses  in  two  pathways  (Fig.  12D;  SI:  posterior  horizontal  layer  ll/lll,  -10±3.7%,  p=0.024,  n=5;  S3: 
posterior  horizontal  layer  V,  -5.6±4%,  p=0.01,  n=6),  facilitated  in  one  pathway  (S2:  anterior  horizontal 
layer  ll/lll,  9.7±2%,  p=0.02,  n=4),  and  had  no  significant  effect  on  the  fourth  pathway  (S4:  vertical  layer 
ll/lll,  p=0.14,  n=7).  Tangential  negative  fields  (-8  V/m,  oriented  anterior  to  posterior)  facilitated  responses 
in  all  horizontal  pathways  (SI:  posterior  horizontal  layer  ll/lll,  11.4±2.5%,  p=0.007,  n=5;  S2:  anterior  layer 
ll/lll,  -10.3±3%,  p=0.009,  n=4;  S3:  posterior  horizontal  layer  V,  4.3±6.9%,  p=0.024,  n=6),  with  no 
significant  modulation  of  synaptic  efficacy  in  the  vertical  pathway  (S4:  vertical  layer  ll/lll,  p=0.13,  n=7). 
Across  modulated  afferent  pathways,  our  results  are  consistent  with  tangential  EFs  that  produce  synaptic 
terminal  hyperpolarization/depolarization  leading  to  facilitation/inhibition  of  synaptic  efficacy  (compare 
Figures  11 C  and  12D).  Though  individual  pathways  were  generally  as  sensitive  to  tangential  fields  as 
radial  fields,  tangential  EFs  had  no  average  effect  on  synaptic  efficacy  averaging  across  pathways  by 
virtue  of  anterior  and  posterior  directed  afferents  (Fig.  12D).  Afferent  axon  terminal  polarization  by 
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subthreshold  DC  fields  was  next  explored  using  compartment  models  to  understand  morphological 
factors  influencing  terminal  polarization. 


Polarization  of  axon  terminals  in  compartment  neuron  models 

Synaptic  terminal  polarization  by  a  uniform  extracellular  EF  leads  to  changes  in  synaptic  efficacy, 
consistent  with  observations  in  cortical  brain  slices.  Therefore,  we  estimated  the  polarization  of  axonal 
terminals  to  weak  DC  stimulation  using  compartment  and  analytical  models.  The  sensitivity  of  a  given 
compartment  to  a  sub-threshold  field  can  be  expressed  as  the  coupling-constant  (or  polarization-length) 
in  terms  of  membrane  polarization  per  unit  EF  (mV  per  V/m)  (Bikson  et  al.,  2004;  Radman  et  al.,  2009a). 

The  relative  terminal  polarization  was  calculated  numerically  for  fields  oriented  parallel  (radial  component 
of  EF,  Fig.  1C2)  or  perpendicular  (tangential  component  of  EF,  Fig.  1C3)  to  the  somatodendritic  axis,  as 
well  as  for  fields  oriented  in  the  optimal  angle  that  leads  to  maximal  polarization  at  each  terminal  (Fig. 
13A,  B).  Layer  V  neurons  with  passive  properties  exposed  to  a  strong  uniform  electric  field  produces  a 
bimodal  polarization  profile  (Fig.  9C2,  false  color  indicates  relative  depolarization  (red)  and 
hyperpolarization  (blue)).  Flowever,  layer  ll/lll  neurons  demonstrate  a  complex  polarization  profile  with 
distal  axonal  branches  approaching  maximal  polarization  at  the  terminals  independently  of  the 
polarization  along  the  neuronal  axis  (Fig.  9C2-3).  Cortical  axons  are  not  straight  even  over  microscopic 
scales  but  we  considered  the  length  of  each  terminating  branch  as  the  distance  from  the  terminal  to  the 
bend  (Fig.  13A).  The  final  branch  space  constant  varied  with  fibre  radius.  We  considered  relative  terminal 
polarization  {Vj  normalized  by  the  axonal  length  constant  and  by  EF)  to  indicate  sensitivity  compared  to 

the  semi-infinite  axon  case  (maximal  polarization=£'A).  For  radial  and  tangential  EFs,  numerical 
simulations  of  cortical  neuron  polarization  indicated  that  terminal  polarization  was  a  complex  function  of 
neuronal  morphology  and  angle  relative  to  the  EF.  Though  varying  within  ±1  (i.e.  ±EX),  there  was  no 
apparent  monotonic  relationship  between  terminal  coupling  constant  and  branch  length  (Fig.  13A,  B1-2). 
Flowever,  further  consideration  of  branch  angle  reveals  that,  specifically  for  long  branches,  relative 
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terminal  polarization  approaches  ±1  for  branches  oriented  parallel  to  the  EF  and  ~0  for  branches 
perpendicular  to  the  EF: 


Equation  3. 


El  if  0  =  0° 
y,(0,f»A)  =  Jo  if  0  =  90° 
-EX  if  0  =  180° 


Moreover,  when  considering  the  optimal  polarization  angle  the  relative  axon  terminal  polarization 
asymptotically  approaches  magnitude  1  with  increasing  branch  length  (Fig.  13B3). 


The  analytical  solution  for  a  bent  axon  (Eq.  1)  was  derived  from  the  cable  equation  to  describe  the 
relationship  between  axon  morphology  and  terminal  coupling  constant.  For  branched  axons  (Fig.  13C1) 
the  terminal  polarization  is  the  sum  of:  a)  The  polarization  at  the  branch  point  (Vg)  weighted  by  the 

inverse  of  the  final  segment’s  electrotonic  length  and,  b)  EX ,  considering  the  component  of  EF  along  the 
final  segment,  weighted  by  the  final  segment’s  electronic  length.  For  short  branches  the  terminal 
membrane  potential  is  coupled  with  the  membrane  potential  at  the  bend  (Vg)  and  therefore  depends  on 
overall  neuronal  morphology  and  where  the  segment  connects  to  the  overall  structure  (Fig.  13B1,  B2, 
note  variability  for  short  branch  lengths).  For  long  branches,  where  the  terminal  membrane 

potential  becomes  independent  of  the  branch  point  and  approaches  EXcosQ .  Equation  1  (using 
numerically  determined  branch  point  voltage)  was  verified  for  numerically  calculated  terminal  polarization 
(Fig.  13D).  For  the  case  of  final  segments  oriented  parallel  to  the  EF,  equation  2  predicts  an  asymptotic 
approach  to  EX  for  final  branch  lengths  >3X,  consistent  with  numerical  simulations  (Fig.  13B3). 
Equation  2  thus  explains  how  overall  neuron  geometry  vs.  axon  terminal  morphology  effects  resulting 
terminal  polarization,  and  how  maximum  expected  polarization  at  a  terminal  would  be  EX  for  a 
sufficiently  long  final  branch  oriented  along  the  field  (Fig.  13B3,  C3).  Importantly,  given  the  density  and 
diversity  of  afferent  axons  in  the  cortex  (some  with  long  optimally-oriented  final  branches)  it  is  reasonable 
to  predict  some  axon  terminals  reach  the  maximal  polarization  EA  for  any  given  EF  orientation. 
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Figure  4.  Modulation  of  pathway-specific  synaptic  efficacy  by  radial  and  tangential  DC  fields. 
Application  of  DC  currents  in  cortical  slice  demonstrates  that  tangential  current  are  as  effective  in 
modulating  pathway-specific  synaptic  efficacy  as  radial  currents,  though  pathway-average  effects 
result  only  for  radial  electric  fields.  A,  Input-output  curve  of  fEPSP  response  amplitude  and  peak 
latency  in  the  horizontal  layer  V  pathway.  Horizontal  grey  bars  indicates  25*^  and  75"^  percentile  of 
fEPSP  peak  latency.  B,  Relative  fEPSP  amplitude  in  the  vertical  layer  V  to  ll/lll  pathway  at 
different  radially  oriented  electric  field  intensities  (correlation  coefficient  of  linear  fit=0.96).  The 
fEPSP  waveform  inset  shows  a  characteristic  change  in  fEPSP  amplitude  with  positive  (+8  V/m, 
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red)  and  negative  radial  fields  (-8  V/m,  blue)  from  control  (no  field,  black).  C,  fEPSP  responses 
are  significantly  (P  <  0.05,  *)  facilitated  with  +8  V/m  fields  (left)  and  reduced  with  -8  V/m  (right)  in 
three  pathways.  In  each  pathway,  individual  slice  averages  are  indicated  with  colored  circles. 
Grouped  average  fEPSP  amplitudes  across  all  synaptic  pathways  indicate  a  7%  polarity-specific 
modulation  of  synaptic  efficacy  with  8  V/m  radial  fields.  Circles  in  the  grouped  average  represent 
the  across  slice  means  of  distinct  pathways  (blue,  red,  green,  and  yellow  are  S1,  S2,  S3,  and  S4 
pathways,  respectively).  D,  fEPSP  amplitude  was  significantly  modulated  by  tangential  electric 
fields  in  all  three  horizontal  pathways  but  with  direction  sensitivity  and  not  in  the  vertical  pathway, 
all  consistent  with  terminal  polarization.  Although  tangential  fields  affected  individual  pathways, 
grouped  average  of  fEPSP  amplitudes  across  all  pathways  was  not  significant. 


The  relevance  of  the  idealized  straight  semi-infinite  axon  (Fig.  1 3C2),  which  can  be  described  by  a  simple 
analytical  solution  (V)  =  EXcosO  or  EX  for  fields  oriented  along  the  branch),  to  realistic  tortuous  cortical 
axon  morphologies  can  now  be  explained  with  the  relationship  we  have  derived  (see  above).  This 
approximation  is  relevant  for  the  cases  of  the  branched  axon  where  the  final  semi-straight  segments’ 
electrotonic  length  is  sufficiently  long  (>  3A);  these  are  also  the  most  sensitive  terminals  to  polarization 
by  EFs.  Thus  the  membrane  length  constant  X  reflects  the  maximal  coupling  constant  of  afferent  axonal 
terminals. 


A  Maximum  Terminal  Polarization 


6  Optimal  polarization  angle  (°) 

Vt  Maximum  terminal  polarization 
at  optimal  angle  (mV  per  mV/mm) 
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Figure  5.  Terminal  polarization  by  uniform  DC  electric  fields  using  neuron  compartment  model  and 
analytical/hybrid  approximations.  Maximum  terminal  polarization  (V))  depends  on  the  length  {£) 
of  the  last  axonal  branch  and  becomes  uncoupled  from  the  bend  point  at  distant  terminals  ( 


26 


DISTRIBUTION  A:  Distribution  approved  for  pubiic  reiease. 


I  >  3A ),  however  for  short  branches  the  terminal  membrane  potential  is  coupled  with  the 
membrane  potential  at  the  bend  {Vq).  In  all  cases,  numerical  simulations  applied  1  V/m  electric 
fields.  A,  For  a  typical  cortical  pyramidal  neuron,  the  maximum  terminal  polarization  (VJ  is  a 
plotted  with  the  corresponding  optimal  polarization  angle  {6)  of  the  branch  relative  to  the  electric 
field  and  the  length  (^)  of  terminating  axon  branch.  B1-2,  Relative  terminal  polarization  (Vt 
normalized  by  the  axonal  length  constant  and  by  the  electric  field)  as  a  function  of  branch 
electrotonic  length  and  angle  (circle  color).  B3,  Considering  the  optimal  polarization  angle  the 
relative  polarization  asymptotically  approaches  magnitude  1  with  branch  length  (equivalently, 
terminal  polarization  reaches  the  maximal  polarization  EX  for  increasing  lengths).  C1-2, 
Schematic  of  a  branched  and  straight  axon  in  a  uniform  electric  field  with  analytical  solutions  (see 
Methods).  The  straight  axon  is  a  special  case  of  the  branched  axon  with  infinite  final  branch 
length.  For  long  branches,  where  ^  »  A  ,  the  terminal  membrane  potential  becomes  independent 
of  the  branch  point  and  approaches  EXcosO .  C3,  The  branched  axon  model  approaches 
maximal  terminal  polarization  EX  for  i>  3X .  D,  Error  of  approximations  (analytical  vs. 
numerical  estimates)  for  branched  (blue)  and  straight  (red)  axons. 


Discussion 

The  application  of  tDCS  in  both  clinical  and  cognitive-neuroscience  research  has  been  encouraged  by  the 
simplicity  of  the  technique  (two  electrodes  and  a  battery  powered  stimulator)  and  the  perception  that 
tDCS  protocols  can  be  designed  by  placing  the  anode/cathode  over  the  cortex  to  “excite/inhibit”.  The 
mechanistic  question  of  how  acute  DCS  modulates  excitability  and  ongoing  neuronal  activity  depends  on 
the  affected  cellular  targets. 


We  characterized  DCS  generated  electric  fields,  how  they  polarize  cellular  compartments,  and  in  turn 
modulate  synaptic  efficacy;  our  approach  applied  the  quasi-uniform  assumption  (Bikson  et  al.,  2012b)  to 
link  finite  element  models  of  cortical  current  flow,  brain  slices,  and  neuron  compartment  models  (Fig.  9). 
On  the  basis  of  our  findings,  the  over-arching  framework  we  propose  is:  During  tDCS  current  crosses  the 
cortex  with  an  electric  field  vector  containing  dominantly  tangential  but  also  radial  components  (relative  to 
the  cortical  surface);  whereas  the  radial  component  will  facilitate/inhibit  synaptic  efficacy  through  somatic 
depolarization/hyperpolarization,  the  tangential  component  will  concurrently  facilitate/inhibit  synaptic 
efficacy,  in  a  pathway  specific  manner,  through  terminal  hyperpolarization/depolarization.  The  rationale 
for  this  framework  in  context  of  prior  mechanistic  studies  is  considered. 
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Animal  studies  of  DCS  on  spontaneous  and  evoked  cortical  activity  (Creutzfeldt  et  al.,  1962;  Purpura  and 
McMurtry,  1965a;  Jefferys,  1981;  Chan  and  Nicholson,  1986;  Chan  et  al.,  1988;  Gluckman  et  al.,  1996; 
Fritsch  et  al.,  2010;  Reato  et  al.,  2010)  demonstrate  that  typically  surface-anode  stimulation  increases 
activity  while  surface-cathode  decreases  activity,  consistent  with  somatic  membrane  polarization 
(Radman  et  al.,  2009a).  Similar  principles  apply  in  the  hippocampus  after  considering  the  inverted 
pyramidal  neuron  morphology  (Gluckman  et  al.,  1996;  Bikson  et  al.,  2004;  Kabakov  et  al.,  2012).  In 
classical  animal  studies,  by  virtue  of  the  stimulation  technique  (e.g.  invasive  electrode)  and/or  cortical 
anatomy  (e.g.  reduced  cortical  folding),  electric  fields  normal  to  the  cortical  surface  were  typically 
generated  (Bindman  et  al.,  1964;  Purpura  and  McMurtry,  1965a).  In  this  vein,  we  confirm  that  purely 
inward  currents  (normal  to  the  cortical  surface)  during  anodal/cathodal  DC  stimulation  in  rat  cortical  slices 
results  in  facilitation/inhibition  of  synaptic  efficacy  (Fig.  11C,  12C);  notably  even  afferents  to  the  apical 
dendrites  were  modulated  based  on  somatic  polarization  (c.f.  (Bikson  et  al.,  2004)).  Specifically, 
polarization  along  the  somatodendritic  axis  with  radial  fields  modulated  cortical  synaptic  efficacy  of 
orthodromically  activated  horizontal  and  vertical  afferents,  consistent  with  postsynaptic  soma  polarization 
(Fig.  lie,  12C)  and  in  agreement  with  results  in  hippocampal  slices  (Jefferys,  1981;  Bikson  et  al.,  2004; 
Kabakov  et  al.,  2012). 

Importantly,  in  the  human  case  with  cortical  current  flow  with  both  a  radial  and  tangential  component;  one 
cannot  ignore  the  tangential  component  and  a  priori  assume  modulation  simply  based  on  somatic 
polarization. 

Our  modeling  results  (Fig.  10)  build  on  previous  efforts  (Datta  et  al.,  2009;  Salvador  et  al.,  2010)  showing 
that  tDCS  generates  both  radial  and  tangential  fields.  But  surprisingly,  quantification  on  macro  and  micro¬ 
scales  demonstrates  tangential  fields  dominate  radial  fields  (4-1 2x)  even  under  the  electrodes  (the 
nominal  stimulation  target).  Most  of  the  cortex  will  have  both  a  radial  and  a  /arger  tangential  component  to 
the  current  flow  (Fig.  10B,  C).  Tangentially  oriented  fields  do  not  significantly  polarize  pyramidal  neuron 
somas  (Ranck,  1975;  Bikson  et  al.,  2006;  Radman  et  al.,  2009a).  Flowever,  tangential  EFs  are  expected 
to  polarize  neuronal  elements  extending  parallel  to  the  cortical  surface,  notably  axons  and  their  synaptic 
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terminals,  including  corticocortical  afferents  (Basser  and  Roth,  2000;  Bikson  et  al.,  2004;  Bikson  et  al., 
2006). 

Generally,  the  sensitivity  of  axon  terminals  to  EFs  and  the  role  of  terminal  membrane  potential  in  synaptic 
function  have  been  recognized  for  decades  (Purpura  and  McMurtry,  1965a;  Hause,  1975;  Bikson  et  al., 
2004;  Bikson  et  al.,  2006).  Though  the  polarization  of  axons  and  their  terminals  in  the  cortex  is  a  complex 
function  of  neuronal  morphology  relative  to  the  EF  (Fig.  9C,  13A),  we  show  sufficiently  long  branch 
terminals  approach  the  polarization  expected  for  semi-infinite  axons:  E’Acos^  (Fig.  13B3,  C3).  If  one 
considers  the  diversity  and  density  of  axonal  afferents,  it  is  reasonable  to  assume  that  some  axons  will 
have  morphology  and  direction  leading  to  optimal  {EX)  terminal  polarization  -  for  any  direction  EF.  X  is 
thus  the  effective  coupling  constant  for  the  most  sensitive  terminals;  -0.42-0. 55  mm  for  cortical  pyramidal 
axons,  (Sasaki  et  al.,  2012).  Cortical  and  hippocampal  pyramidal  somas  have  a  typical  coupling  constant 
-0.12-0. 24  mm  (Jefferys,  1981;  Bikson  et  al.,  2004;  Deans  et  al.,  2007;  Radman  et  al.,  2009a;  Frohlich 
and  McCormick,  2010).  This  directionless  terminal  coupling  constant  is  independent  of  and  2-3  times 
more  sensitive  than  cortical  pyramid  somas  under  optimally-oriented  radial  EFs  (Flause,  1975;  Radman  et 
al.,  2009a).  Thus,  even  under  purely  radial  fields,  afferent  axon  terminals  will  polarize  more  than  pyramid 
somas.  Moreover,  during  tDCS  EFs  are  dominantly  tangential  (Fig  10),  which  might  indicate  a  further 
reduction  of  pyramid  soma  polarization,  but  not  afferent  axon  terminals. 

The  direction  of  terminal  polarization  will  depend  on  the  orientation  (morphology)  of  the  specific  afferent 
pathway  relative  to  the  EF  (Fig.  11C,  12D).  This  is  especially  relevant  for  the  cortex,  where  different  types 
of  excitatory  and  inhibitory  cells  form  dense  anatomical  and  functional  synaptic  circuits  that  establish 
routes  of  information  transfer.  As  a  result  of  the  symmetrical  orientation  of  the  local  synaptic  pathways  in 
the  primary  motor  cortex,  we  demonstrate  tangential  EFs  have  no  average  effect  on  synaptic  efficacy 
despite  pathway  specific  effects  (Fig.  12D)  -  importantly,  given  distinct  pathway  specific  functions,  a 
negligible  average  effect  does  nof  equate  to  no  net  effect  on  information  processing.  Concomitant  radial 
and  tangential  field  effects  in  vivo  should  be  addressed  in  future  studies.  Our  results  in  cortical  brain 
slices  are  consistent  with  tangential  EFs  that  produce  synaptic  terminal  hyperpolarization/depolarization 
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leading  to  facilitation/inhibition  of  synaptic  efficacy  (Fig.  11 C,  12D).  Similar  directional  sensitivity  to  electric 
fields  oriented  along  afferent  axons  is  reported  in  the  hippocampus  (Bikson  et  al.,  2004;  Kabakov  et  al., 
2012).  These  results  are  in  agreement  with  previous  findings  on  the  role  of  synaptic  terminal  polarization 
on  synaptic  strength,  as  presynaptic  hyperpolarization  increases  presynaptic  action  potential  size  and 
modulates  Ca^'"  current  activation  and  driving  force  (Hubbard  and  Willis,  1962b;  Takeuchi  and  Takeuchi, 
1962;  Miledi  and  Slater,  1966).  Accordingly  release  is  enhanced  by  terminal  hyperpolarization  and 
decreased  by  terminal  depolarization  (Bullock  and  Hagiwara,  1957;  Hubbard  and  Willis,  1962a;  Takeuchi 
and  Takeuchi,  1962;  Miledi  and  Slater,  1966;  Hubbard  and  Willis,  1968;  Dudel,  1971). 


Conclusion 

We  developed  an  overarching  framework  for  modulation  of  synaptic  efficacy  by  acute  DCS  based  on  a 
simple  qualitative  model  where  the  changes  in  synaptic  efficacy  are  determined  by  somatic  and  afferent 
axon  polarization  (Fig.  11C).  The  results  we  observed  across  four  pathways  for  all  field  orientation  and 
polarities  were  consistent  with  this  model  (Fig  12C,  D).  Synaptic  terminals,  therefore,  should  also  be 
considered  as  an  additional  target  of  modulation  during  DC  stimulation,  provided  the  evidence  for  a 
presynaptic  contribution  to  DC  stimulation  effect  (Kabakov  et  al.,  2012;  Marquez-Ruiz  et  al.,  2012)  and 
especially  given  the  role  of  terminals  in  information  processing  and  plasticity  (Malenka  and  Nicoll,  1999). 
Indeed,  in  this  last  sense  the  role  of  somatic  vs.  terminal  effects  can  be  considered  analogous  to  the  long¬ 
standing  debate  over  pre-  and  post-synaptic  locus  of  long-term  potentiation  and  depression  of  synaptic 
efficacy.  Both  will  likely  prove  important  since  during  tDCS  cortical  regions  are  exposed  to  EFs  with  both 
radial  and  tangential  components,  such  that  the  somatic  effects  produced  by  the  radial  component  would 
be  modulated  by  concurrent  afferent-specific  terminal  effects  in  an  antagonistic  or  complementary 
fashion. 

The  cellular  locus  of  acute  DCS-induced  excitability  changes  (e.g.  somas,  dendrites,  synaptic  terminals) 
is  fundamental  to  understanding  the  mechanisms  of  tDCS  and  rationalizing  stimulation  strategies.  While 
this  study  establishes  the  principle  of  axon  terminal  effects,  alongside  somatic  effects,  as  a  plausible 
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target  of  modulation  by  DC  stimulation,  it  remains  to  be  established  where  the  functional  outcomes  of 
DCS  (tDCS)  originate.  For  example,  this  study  addressed  only  the  acute  synaptic  and  non-synaptic 
effects  of  short  duration  DC  stimulation,  thus  we  did  not  expect  long-term  plastic  changes.  Additionally, 
interneuronal/inhibitory  effects  on  synaptic  function  were  excluded  in  the  present  study,  though  actions  on 
terminals  indicate  a  potential  for  modulation  if  inhibitory  function,  even  if  interneuron  somata  are  not 
significant  polarized  (Radman  2009).  Moreover,  brain  slices  are  relatively  quiescent  compared  to  the  in 
situ  case,  with  no  oscillatory  activity  or  stimulation-matched  activation  -  these  factors  can  be  further 
systematically  addressed  using  brain  slice  techniques  (Nagarajan  et  al.,  1993;  Frohlich  and  McCormick, 
2010;  Reato  et  al.,  2010).  Notwithstanding  these  important  features,  understanding  the  direct  modulation 
of  cortical  excitatory  transmission  by  acute  DCS  is  a  necessary  building  block  toward  a  broader 
mechanistic  understanding. 
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Section  2:  Diffuse  and  sustained  electrical  stimulation  amplifies  synaptic 

cooperativity 


Abstract 

Transcranial  direct  current  stimulation  (tDCS)  produces  weak  and  diffuse  current  flow  in  the  brain  with 
effects  that  are  state-dependent  and  outlast  stimulation  when  current  is  sustained.  A  mechanistic 
framework  for  tDCS  should  capture  these  spatiotemporal  features.  It  remains  unclear  how  electrical 
stimulation  affects  ongoing  synaptic  dynamics  and  how  modulation  of  afferent  inputs  changes  synaptic 
activity  at  the  target  brain  region.  We  tested  the  effect  of  acute  direct  current  stimulation  (DCS;  10-20  V/m 
for  3-5  s)  on  constant  rate  (5-40  Hz)  and  Poisson-distributed  (4  Hz  mean)  trains  of  synaptic  activity  during 
periods  of  synaptic  adaptation  or  post-adaptation.  Sustained  synaptic  activity  during  DCS  was  polarity- 
specifically  modulated  across  tested  input  frequencies.  Synaptic  depression  attenuates  the  sensitivity  to 
DCS  from  1.1%  per  V/m  to  0.55%.  Sustained  DCS  thus  facilitates  cumulative  neuromodulation, 
potentially  reversing  endogenous  synaptic  depression.  We  establish  these  effects  are  mediated  by 
afferent  axon  fiber  polarization,  which  boosts  cooperativity  between  synaptic  inputs  when  the  electric  field 
is  oriented  along  the  fibers.  This  potentially  extends  the  locus  of  neuromodulation  from  the  nominal  brain 
target  to  afferent  brain  regions.  Based  on  these  results  we  hypothesized  that  the  polarization  of  afferent 
neurons  in  upstream  brain  regions  may  modulate  activity  in  the  target  brain  region  during  tDCS.  A 
multiscale  model  of  electrical  stimulation  including  a  finite  element  model  of  brain  current  flow,  numerical 
simulations  of  neuronal  activity,  and  a  statistical  theory  of  coincident  activity  predicts  the  diffuse  and  weak 
profile  of  current  flow  can  be  advantageous. 


Introduction 

Research  into  the  cognitive  and  behavioral  consequences  of  transcranial  Direct  Current  Stimulation 
(tDCS)  has  outpaced  development  of  cellular  models  that  can  explain  the  diversity  of  applications, 
especially  in  the  context  of  the  unique  spatiotemporal  features  of  tDCS:  sustained  flow  of  weak  direct 
current  across  large  areas  of  the  brain.  Computational  models  of  brain  current  flow  during  tDCS  predict 
low  intensity  and  diffuse  current,  which  weakly  polarizes  neurons  across  the  brain.  Clinical 
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neurophysiology  and  imaging  has  demonstrated  that  weak  and  diffuse  current  flow  produces  online 
changes  in  cortical  excitability  that  is  sustained  for  the  duration  of  stimulation.  It  remains  unclear  how 
diffuse  and  low  intensity  currents  can  influence  connectivity  within  a  local  circuit  and  between  brain 
regions. 

Low-intensity  direct  current  stimulation  (DCS)  produces  a  small  change  in  membrane  potential  (<  1  mV, 
(Radman  et  al.,  2009a)),  which  is  not  sufficient  to  induce  spiking  activity  in  resting  cortical  pyramidal  cells. 
DCS,  while  weak,  modulates  synaptic  efficacy  at  rest  as  measured  by  single  evoked  responses  in 
humans  ((Nitsche  and  Paulus,  2000;  Priori,  2003))  as  well  as  in  animals  ((Bikson  et  al.,  2004;  Kabakov  et 
al.,  2012;  Marquez-Ruiz  et  al.,  2012;  Rahman  et  al.,  2013)).  tDCS  is  typically  applied  adjunct  to  a  task 
(corresponding  to  sustained  neuronal  activity)  in  clinical  and  performance  enhancement  applications.  In 
these  applications  both  the  DCS  and  the  endogenous  activity  are  sustained.  However,  the  role  of  DCS 
acutely  on  sustained  synaptic  activity  is  unclear.  In  what  ways  might  the  dynamics  of  endogenous  activity 
enhance  or  dull  the  effects  of  sustained  DCS? 

Acute  synaptic  modulation  by  DCS  has  been  linked  to  both  polarization  along  the  somato-dendritic  axis 
(Jefferys,  1981;  Bikson  et  al.,  2004;  Rahman  et  al.,  2013)  and  afferent  axon  terminals  (Hause,  1975; 
Kabakov  et  al.,  2012;  Boizoni  et  al.,  2013a;  Rahman  et  al.,  2013;  Baczyk  and  Jankowska,  2014). 
Neuronal  polarization  does  not  have  to  be  limited  to  just  the  target  brain  region  during  tDCS;  activity  in  the 
target  brain  region  can  be  modulated  by  the  simultaneous  polarization  of  neurons  and  their  processes  in 
afferent  brain  regions.  Specifically,  synaptic  function  in  the  target  region  will  be  determined  by  both  local 
synaptic  efficacy  and  afferent  activity  in  terms  of  number  of  active  axons  and  their  firing  rate  -  together 
these  set  the  connectivity  between  the  afferent  region  and  the  target.  We  propose,  sustained  DCS 
coupled  with  endogenous  activity  may  modulate  neuronal  response  in  the  target  brain  region  through 
changes  in  afferent  activity  as  a  direct  consequence  of  the  diffuse  current  flow. 

We  quantified  changes  in  synaptic  efficacy  in  rat  brain  slices  during  constant  and  Poisson  rate  afferent 
synaptic  activity.  We  assessed  for  the  first  time  DCS  modulation  of  synaptic  efficacy  during  adaptation  to 
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ongoing  activity  and  post-adaptation.  Our  experiments  identify  a  role  for  afferent  axon  polarization  in 
driving  ongoing  activity.  Experimental  results  are  integrated  into  a  model  of  connectivity  between  brain 
regions  during  tDCS  which  predicts  that  diffuse  and  sustained  current  flow  increases  the  probability  of 
coincident  activity  and  endogenous  plasticity  between  brain  regions.  Our  findings  provide  a  substrate  for 
how  weak  and  diffuse  DCS  can  amplify  synaptic  activity,  modulate  synaptic  dynamics,  and  boost  synaptic 
efficacy  through  cooperativity  of  afferent  inputs. 


Methods 


Electrophysiology 

All  animal  experiments  were  carried  out  in  accordance  with  guidelines  and  protocols  approved  by  the 
Institutional  Animal  Care  and  Use  Committee  (lACUC)  at  The  City  College  of  New  York,  CUNY  (Protocol 
No:  846.3). 

Brain  slices  including  a  part  of  the  primary  motor  cortex  (M1)  were  prepared  from  male  young  adult  Wistar 
rats  aged  3  to  6  weeks  old,  which  were  deeply  anesthetized  with  ketamine  (7.4  mg/kg)  and  xylazine  (0.7 
mg/kg)  applied  intraperitoneally  (i.p.)  and  sacrificed  by  cervical  dislocation.  The  brain  was  removed  and 
immersed  in  chilled  (2-6  °C)  artificial  cerebrospinal  fluid  (ACSF)  containing  (in  mM):  126  NaCI,  4.4  KCI, 
1.25  NaH2P04,  2  MgS04,  2  CaCI,  24  NaHC03  and  10  D-glucose,  bubbled  with  a  mixture  of  95%  02-5% 
C02).  A  slightly  modified  ACSF  composition  was  used  for  long-term  plasticity  experiments  (1  mM 
MgS04).  Coronal  slices  (400  pm  thick)  were  cut  using  a  vibrating  microtome  and  transferred  to  a  holding 
chamber  for  at  least  30  min  in  ambient  temperature.  Slices  were  then  transferred  to  a  fluid-gas  interface 
chamber  perfused  with  warmed  ACSF  (30.0±0.5  °C)  at  1.9  mL/min.  The  humidified  atmosphere  over  the 
slices  was  saturated  with  a  mixture  of  95%  02-5%  C02.  Recordings  started  1-3  h  after  dissection. 

We  stimulated  the  excitatory  vertical  pathway  from  L5A^2/3  and  the  horizontal  pathway  within  the 
superficial  layers  (L2/3).  A  recording  electrode  was  placed  250-400  pm  below  the  pial  surface  in  L2/3  and 
a  bipolar  platinum/stainless  steel  stimulating  electrode  was  positioned  -500  pm  vertically  below  it  at  the 
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boundary  of  L3  and  L5.  Either  15  (in  adaptation  experiments)  or  200  (in  post-adaptation  experiments) 
constant-current  pulses  (0.2  ms,  10-150  pA)  were  delivered  at  5,  10,  20,  or  40  Hz  for  constant  train 
experiments  with  an  interval  of  60  s  (15  pulse  trains)  or  4  min  (200  pulse  trains)  between  trains  to  allow 
for  recovery  from  synaptic  depression.  Stimulus  intensity  was  adjusted  to  produce  responses  half  of  the 
maximum  amplitude  that  could  be  evoked.  In  a  subset  of  experiments  to  evaluate  the  effect  of  DCS  on  a 
natural  train  of  presynaptic  action  potentials,  81  pulses  were  delivered  with  Poisson  distributed  interspike 
intervals  (4  Hz  mean,  20  sec). 

Uniform  extracellular  EFs  (±10  and  20  V/m)  were  generated  by  passing  constant  current  (D/A  driven 
analog  follower;  A-M  Systems,  WA,  USA)  between  two  large  Ag-AgCI  wires  positioned  in  the  bath  across 
the  slice  starting  0.5  s  before  a  train  of  constant  or  Poisson-distributed  stimuli.  In  post-adaptation 
experiments,  the  field  was  applied  after  100  pulses  for  5  sec  that  consisted  of  50  pulses  during-DCS 
followed  by  another  50  pulses  post-DCS  for  a  total  of  200  pulses  at  a  constant  frequency  of  1 0  Hz.  For  all 
experiments,  we  used  a  two-tailed  t-test  and  a  bootstrap  hypothesis  test  to  test  for  significance. 


Modeling  synaptic  dynamics  and  recruitment 

The  neuronal  population  response  to  excitatory  presynaptic  drive  can  be  modeled  as  the  averaged 
voltage  response  7(t)  (Richardson  et  al.,  2005).  A  train  of  presynaptic  spikes  arriving  down  input  fiber  n 
over  a  large  population  of  input  fibers  at  time  t  evokes  excitatory  postsynaptic  potentials  modeled  as 
alpha  synapses  a(t). 


Equation  4.  7(t)  =  “  ^k) 

The  set  of  amplitudes  {4^.},  reflecting  the  field  potential  amplitudes,  are  modeled  using  a 
phenomenological  description  of  dynamic  synapses  with  short-term  plasticity  (Tsodyks  and  Markram, 
1997;  Mongillo  et  al.,  2008).  A  simple  quantitative  model  of  postsynaptic  response  evoked  by  multiple 
synapses  can  be  described  as  a  product  of  a  constant  Ag  and  two  variables  F  and  D: 
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Equation  5.  A{t)  =  Ag  ■  f  (t)  ■  D(t) 


dF{t') 

dt 


Fq  -  F(t) 

Tf 


t  —  tj^p  ^  F  {F)  —^F{t  )  +  A 


dD(t)  l-D(t) 
di  T  p) 


Synaptic  efficacy  is  modulated  by  the  amount  of  available  resources  D  and  the  utilization  parameter  F 
that  defines  the  fraction  of  resources  used  by  each  spike.  The  utilization  parameter  may  reflect  residual 
calcium  level  in  the  presynaptic  terminal.  Upon  a  spike,  an  amount  FD  of  the  available  resources  is  used 
to  produce  the  postsynaptic  current,  thus  reducing  D\  physiologically  the  process  mimics  neurotransmitter 
depletion.  The  spike  also  increases  f,  mimicking  presynaptic  calcium  influx.  Experimentally,  the  field 
potential  amplitude  reflects  the  population  averaged  synaptic  current  and  is  well  approximated  by  this 
model. 


Between  spikes,  F  and  D  exponentially  recover  to  their  baseline  levels  (D  =  1  and  F  =  F^,  F  <  1)  with 
time  constants  Xp  (facilitating)  and  (depressing).  The  phenomenological  model  reproduces  behavior  of 
both  facilitating  {xp  >  Xp,)  and  depressing  (xp,  >  Xp)  cortical  synapses. 

We  simulated  a  single  postsynaptic  integrate-and-fire  neuron  receiving  correlated  excitatory  input  signals. 
The  cell  received  input  through  at  most  1000  synapses  during  ongoing  synaptic  activity  at  a  constant  rate 
(10  Hz).  To  mimic  glutamatergic  transmission,  synapses  were  conductance-based  with  reversal  potential 
=  and  time  constant  Xp^^  =  Sms.  Additionally,  every  presynaptic  spike  leads  to  synaptic 
depression.  Recruitment  was  modeled  as  an  increase  or  decrease  in  the  number  of  active  synapses.  The 
number  of  active  synapses  was  decreased  from  1000  to  500  to  simulate  the  dynamics  during  a  decrease 
in  synaptic  cooperativity  or  increased  from  500  to  1000  to  simulate  dynamics  during  an  increase  in 
synaptic  cooperativity. 
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Statistical  theory  of  coincident  pre-  and  postsynaptic  spikes 

We  present  a  statistical  theory  (Equation  5)  for  increasing  synaptic  strength  through  changes  in  the 
likelihood  of  pre-  or  postsynaptic  firing  by  considering  the  number  of  presynaptic  spikes  {Npre)  ^nd  the 
number  of  postsynaptic  spikes  (iV^ost)  in  discrete  time  bins  (iV^ms).  The  duration  of  each  bin  is  taken  to  be 
10-20  ms  corresponding  to  the  time  constant  of  spike  timing  dependent  plasticity.  The  hypergeometric 
distribution  is  used  to  estimate  the  coincidence  count  probability  P{k\Npre,Npost>Nbins)  of  observing  k 
coincident  pre-  and  postsynaptic  spikes  from  a  given  number  of  bins  and  having  observed 
postsynaptic  spikes  and  presynaptic  spikes. 


Equation  6.  P(fc|iVp^e- 

Results 

A  quantitative  description  of  synaptic  dynamics  during  DCS 

The  effects  of  DCS  on  online  information  processing  in  the  brain  can  only  be  understood  by  considering 
synaptic  dynamics  during  ongoing  activity.  A  fundamental  feature  of  synaptic  transmission  is  adaptation 
to  continuous  afferent  inputs  (Abbott  and  Regehr,  2004),  which  is  quantitatively  described  by  short-term 
plasticity  (Equation  1).  Adaptation  to  synaptic  activity  arising  from  the  depletion  of  neurotransmitters  at  the 
presynaptic  terminal  underlies  important  neural  functions  like  gain  control  (Abbott  et  al.,  1997)  and 
working  memory  (Mongillo  et  al.,  2008).  The  short-term  plasticity  phenomenon  regulates  synaptic  efficacy 
and  is  quantified  as  a  gradual  decrease  (for  depressing  synapses)  or  an  increase  (for  facilitating 
synapses)  in  the  excitatory  postsynaptic  potential  (EPSP,  Figure  14B)  during  continuous  presynaptic 
inputs.  Specifically,  we  quantified  how  DCS  during  presynaptic  activity  leads  to  a  sustained  change  in  the 
EPSP  (Equation  1).  A  conductance-based  neuron  model  with  dynamic  synapses  suggests  that  changes 
in  postsynaptic  membrane  potential  during  DCS  produces  a  sustained  modification  in  EPSP  amplitude 
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and  produces  a  greater  net  depolarization  when  anodal  DCS  is  applied  during  ongoing  synaptic  activity 
(Figure  14C).  To  test  this  prediction,  we  used  an  in  vitro  preparation  where  we  could  control  the  frequency 
of  presynaptic  inputs  and  study  the  effects  of  DCS  on  synaptic  efficacy  through  the  extracellularly- 
recorded  (field)  EPSP  (Figure  15A),  which  we  refer  to  simply  as  the  field  potential  (FP). 
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Figure  6.  The  postsynaptic  voltage  response  during  DCS  and  ongoing  presynaptic  activity  results 
in  sustained  and  cumulative  changes  that  are  regulated  by  synaptic  efficacy,  number  of  active 
inputs,  and  rate  of  presynaptic  activity.  A)  Schematic  representing  the  voltage  output  of  a 
postsynaptic  cell  during  a  train  of  presynaptic  spikes  arriving  at  times  [t^k]  down  input  fibers. 
Synaptic  transmission  at  cortical  synapses  are  not  static,  but  are  dynamically  regulated  by  the 
available  pool  of  releasable  vesicles.  B)  Synaptic  efficacy  during  activity  is  controlled  by  short¬ 
term  plasticity  based  on  the  initial  vesicle  release  probabilities  (Preiease)  and  vesicle  depletion  and 
recovery  times.  C)  Simulation  of  the  postsynaptic  response  during  DCS  in  a  conductance-based 
model  with  dynamic  synapses.  Anodal  DCS  (red),  modeled  as  a  postsynaptic  depolarization, 
facilitates  EPSP  amplitudes  and  cathodal  DCS  (blue),  modeled  as  a  postsynaptic 
hyperpolarization,  depresses  EPSP  amplitudes.  The  postsynaptic  response  during  DCS  is 
sustained  for  the  duration  of  DCS. 


DCS  has  a  sustained  and  cumulative  effect  on  synaptic  efficacy 

The  effects  of  DCS  on  synaptic  efficacy  were  evaluated  during  ongoing  presynaptic  activity  in  the  rat 
primary  motor  cortex.  Presynaptic  afferent  axons  (L5A^2/3  pathway)  were  stimulated  at  5,  10,  20,  and 
40  FIz  to  simulate  synaptic  activity.  Field  potentials  (FPs)  decayed  in  amplitude  over  the  course  of  a  15 
pulse  train  as  a  result  of  synaptic  depression  (Figure  15B).  Anodal  DCS  increased  (+10  V/m:  13.2±6.9%, 


38 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


+20  V/m:  20.3±9.8%,  n=124,  p<0.05)  and  cathodal  DCS  decreased  (-10  V/m:  -13.3±8.7%,  -20  V/m:  - 
23.3±11.2%,  n=124,  p<0.05)  the  transient  FP  amplitude  when  applied  during  the  train  of  inputs  (Figure 
15C).  Importantly,  we  show  for  the  first  time  that  the  change  in  synaptic  efficacy  is  maintained  for  the 
duration  of  DCS  when  there  is  a  constant  rate  of  synaptic  activity  (Figure  15B,  15C).  The  steady-state  FP 
amplitude  during  DCS  is  significantly  greater  during  anodal  DCS  and  smaller  during  cathodal  DCS, 
compared  to  the  control  case  (no-DCS,  Figure  15B,C). 


Computational  models  suggest  there  is  cumulative  membrane  depolarization  during  ongoing  synaptic 
activity  and  DCS  may  amplify  the  net  depolarization  (Figure  14C).  In  vitro,  DCS  has  a  greater  cumulative 
increase  in  synaptic  efficacy  during  activity.  Anodal  DCS  increased  and  cathodal  DCS  decreased  how 
fast  the  cumulative  synaptic  efficacy  changed  during  synaptic  activity,  relative  to  baseline  (Figure  14F, 
15C).  In  slices  with  weak  synaptic  depression  or  at  low-frequencies  of  synaptic  activity,  anodal  DCS  can 
reverse  synaptic  depression  towards  facilitation  for  the  duration  of  DCS.  This  finding  is  remarkable 
because  it  indicates  DCS  can  qualitatively  change  the  nature  of  gain  control  in  a  given  brain  region, 
turning  depressive  dynamics  into  facilitating  ones. 
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Figure  7.  DCS  in  rat  primary  motor  cortex  results  in  a  sustained  modulation  of  synaptic  efficacy.  A) 
Schematic  of  anodal  (left)  and  cathodal  (right)  DCS  with  current  flow  along  the  somatodendritic 
axis  and  the  effective  membrane  polarization  of  cells  and  axons  in  false-color.  B)  Orthodromic 
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stimulation  of  the  L5-*2/3  pathway  at  a  constant  frequency  (20  Hz  shown)  in  M1  transiently 
facilitates  synaptic  efficacy  during  anodal  DCS  and  depresses  efficacy  during  cathodal  DCS.  The 
DCS  effects  are  sustained  for  the  duration  of  stimulation.  C)  DCS  modulates  synaptic  efficacy 
during  continuous  afferent  inputs  despite  synaptic  depression  arising  from  vesicle  depletion  in  the 
presynaptic  terminals.  D)  The  transient  change  in  field  potential  (FP)  amplitude  during  DCS  is 
correlated  with  the  steady-state  change  in  FP  amplitude.  FP  amplitudes  in  the  constant-rate  train 
of  inputs  during  DCS  are  normalized  by  the  corresponding  n-th  FP  amplitude  during  a  train 
without  DCS.  E)  The  sensitivity  (%  change  in  FP  amplitude  per  V/m)  to  DCS  of  the  steady-state 
FP  is  less  than  the  sensitivity  of  the  transient  FP  (n=124).  Gray  best  fit  lines  are  for  each  of  the 
tested  frequencies  (5,  10,  20,  and  40  Hz)  and  red  is  the  combined  average  across  frequencies. 
There  was  no  significant  difference  in  either  transient  or  steady-state  sensitivity  across 
frequencies.  F)  The  rate  of  cumulative  change  in  FP  amplitude  (the  cumulative  %  change  from 
baseline)  is  greater  during  anodal  DCS  and  less  during  cathodal  DCS  compared  to  control.  Dark 
points  indicate  20  V/m  and  light  points  are  10  V/m.  In  all  cases,  magenta  indicates  anodal  DCS 
and  blue  is  cathodal  DCS. 


DCS  produces  frequency-independent  changes  in  synaptic  efficacy 

A  model  of  dynamic  synaptic  transmission  consisting  of  a  facilitating  (F)  and  a  depressing  (D)  term  with 
first-order  decay  kinetics  (referred  to  as  the  FD  model)  fit  to  the  data  predicted  the  observed  EF  effects. 
The  most  common  form  of  short-term  plasticity  (>82%  of  experiments)  observed  in  the  rat  primary  motor 
cortex  was  short-term  depression  (t^  >  v,  =  0.74  sec,  Tp  =  0.23  sec,  f,,  =  0.18,  A  =  0.41).  DCS  had  no 
effect  on  the  time  constants  for  depression  or  facilitation  (Wilcoxon  signed-rank  test,  p  >  0.05),  suggesting 
DCS  does  not  have  a  direct  effect  on  presynaptic  transmitter  release  kinetics.  Both  in  our  experiments 
and  in  the  FD  model,  ongoing  synaptic  activity  results  in  a  substantial  reduction  in  FP  amplitude  due  to 
synaptic  depression.  The  change  in  synaptic  efficacy  resulting  from  DCS  (and  independent  of  synaptic 
depression)  is  quantified  by  normalizing  each  n-th  FP  amplitude  during  DCS  (FPn.ocs)  by  the 
corresponding  n-th  FP  amplitude  without  DCS  (FPn, control),  within-slice.  The  normalization  reveals  a 
sustained  increase  of  synaptic  efficacy  with  anodal  DCS  and  decrease  with  cathodal  DCS  for  as  long  as 
the  EF  is  on  (Figure  16A).  The  sustained  change  in  synaptic  efficacy  is  frequency-independent  and 
persists  during  a  naturalistic  pattern  of  synaptic  activity  (Figure  17). 
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Figure  8.  Synaptic  efficacy  is  frequency-independently  modulated  during  DCS,  sustained  for  the 
duration  of  stimulation,  and  is  cumulative.  A)  FP  amplitudes  during  20  V/m  EFs  are  normalized  by 
the  corresponding  n-th  FP  amplitude  in  a  control  train,  within  slice.  B)  Anodal  DCS  attenuates 
synaptic  depression.  The  cumulative  percent  change  from  baseline  is  plotted  for  anodal 
(magenta),  cathodal  (blue),  and  control  (black)  conditions.  C)  The  transient  change  in  FP 
amplitude  is  correlated  with  the  steady-state  change  in  FP  amplitude. 
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Figure  9.  DCS  produces  a  sustained  change  in  synaptic  efficacy  during  a  naturalistic  pattern  of 
presynaptic  activity.  Afferent  inputs  reflecting  the  Poisson  distributed  spike  train  pattern  in  vivo 
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(mean  rate  4  Hz)  is  acutely  facilitated  by  anodal  DCS  (magenta)  and  depressed  by  cathodal  DCS 
(blue). 

Ongoing  synaptic  activity  changes  the  synaptic  sensitivity  to  DCS 

There  is  a  marked  decrease  in  the  sensitivity  to  DCS  over  time,  quantified  as  the  %  change  in  FP 
amplitude  from  control  per  V/m,  from  1.17±0.1%  /  V/m  to  0.65±0.05%  /  V/m  during  ongoing  synaptic 
activity  (Figure  1 5E).  This  change  in  sensitivity  is  independent  of  the  rate  of  synaptic  activity  (Figure  1 5E, 
gray  lines  indicate  different  rates  of  afferent  synaptic  inputs;  rate  coefficients:  transient:  0.03±0.06%  /  V/m; 
steady-state:  -0.03±0.07%  /  V/m).  Furthermore,  the  transient  and  steady-state  change  in  synaptic  efficacy 
is  directly  correlated,  across  all  frequencies  (Figure  15D,  Figure  16B).  Synaptic  adaptation,  therefore, 
reduces  the  capacity  for  DCS  to  modulate  synaptic  efficacy. 


DCS  modulates  synaptic  efficacy  post-adaptation 

Modulation  of  synaptic  efficacy  after  the  adaptation  period  was  measured  by  applying  a  long-duration 
train  of  200  stimuli  at  10  Hz  (Figure  18A).  20  V/m  EFs  were  applied  for  5  sec  after  the  first  100  pulses  in 
the  train.  During  stimulation,  synaptic  efficacy  transiently  decreased  with  cathodal  DCS  (-33.5±21%,  n=9, 
p<0.05)  and  increased  with  anodal  DCS  (28.8±26.8%,  n=11,  p<0.05).  The  effect  of  DCS  was  sustained 
for  the  duration  of  stimulation  but  there  is  an  adaptation  to  DCS  during  anodal  stimulation  (Figure  18B). 


DCS  promotes  synaptic  cooperativity  through  recruitment  of  afferent  axons 

The  negative  correlation  between  the  transient  change  in  FP  amplitude  and  the  change  in  FP  amplitude 
during  DCS  suggests  there  is  an  adaptation  to  the  effect  of  DCS  (Figure  18C).  We  hypothesized  the 
adaptation  to  DCS  may  be  attributed  to  recruitment  of  presynaptic  axons  during  DCS.  EFs  along  the 
ascending  afferent  axons  may  be  polarizing  presynaptic  compartments  leading  to  recruitment  of 
additional  afferents  {N^  in  Equation  1)  during  orthodromic  stimulation.  We  tested  this  hypothesis  by 
changing  the  orthodromic  stimulation  intensity  by  >2X  or  by  2x,  thus  decreasing  or  increasing. 
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respectively,  the  number  of  afferent  fibers  generating  the  field  potential.  The  protocol  is  similar  to  the 
post-adaptation  DCS  experiment  but  here  the  EF  is  replaced  with  a  change  in  stimulus  intensity. 

Both  DCS  and  changes  in  orthodromic  stimulus  intensity  result  in  similar  dynamics  of  the  FP  (Figure 
18D).  During  cathodal  DCS  and  the  analogous  Yzx  stimulus  intensity  experiments,  FP  amplitudes  were 
transiently  depressed  and  recovered  from  depression.  During  anodal  DCS  and  the  analogous  2x 
orthodromic  stimulus  intensity  experiments,  FP  amplitudes  were  transiently  facilitated  and  gradually 
recovered  from  facilitation  towards  baseline.  In  both  cathodal  DCS  and  >2X  stimulus  intensity  experiments 
we  observe  a  prominent  rebound  in  FP  amplitude  after  DCS  was  turned  off  and  after  stimulus  intensity 
was  returned  to  the  probing  intensity.  These  results  suggest  DCS  changes  the  presynaptic  axon 
membrane  potential  and  modulates  the  sensitivity  to  firing  action  potentials.  A  model  incorporating 
synaptic  depression  and  recruitment  is  able  to  capture  the  salient  features  of  the  experimental  data. 
Recruitment  is  modeled  as  an  instantaneous  increase  in  the  number  of  active  presynaptic  inputs  for  the 
duration  of  stimulation.  A  DCS-induced  postsynaptic  depolarization  alone  fails  to  capture  the  nonlinear 
effects  of  DCS  on  FP  amplitude. 

DCS  can  change  synaptic  efficacy  even  without  recruitment  of  afferent  fibers  during  orthodromic 
stimulation.  Axons  within  L2/3  (L2/3^2/3)  were  orthodromically  stimulated  while  only  polarizing  the 
postsynaptic  cells.  There  is  a  significant  transient  change  in  FP  amplitude  during  DCS  (-20  V/m: 
3.15±1.18%,  n=5;  +20  V/m:  6.44±4.78%,  n=5).  The  dynamics  of  recruitment  which  would  have  resulted  in 
an  offset  of  the  FP  amplitude  and  a  recovery  from  facilitation  or  depression  towards  baseline  were  not 
observed.  This  result,  along  with  previously  published  work,  shows  postsynaptic  polarization  alone  can 
modulate  synaptic  efficacy.  Flowever,  the  transient  change  in  synaptic  efficacy  is  less  than  when 
activating  the  vertical  L5^2/3  pathway  (6.44%  vs.  28.8%  for  20  V/m)  suggesting  some  of  the  change  in 
synaptic  efficacy  we  observe  during  stimulation  of  the  L5^2/3  pathway  may  result  from  recruitment  of 
afferent  axons. 
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Figure  10.  DCS  post-adaptation  modulates  synaptic  efficacy  and  is  driven  by  recruitment  of  afferent 
axons.  A)  Cathodal  (blue)  and  anodal  DCS  (magenta)  during  a  train  of  synaptic  inputs  at  10  Hz 
modulates  synaptic  efficacy  for  the  duration  of  stimulation  despite  adaptation.  B)  DCS  polarity 
specifically  modulates  the  transient  FP  amplitude  (first  evoked  potential  during  DCS)  relative  to 
pre-DCS.  C)  There  is  a  correlation  between  the  transient  change  in  FP  amplitude  and  the  amount 
of  recovery  from  transient  facilitation  or  depression  during  DCS.  D)  Changing  the  stimulus 
intensity,  instead  of  applying  DCS,  results  in  similar  dynamics  on  FP  amplitude.  Decreasing 
stimulus  intensity  by  a  half  (left)  decreases  FP  amplitudes.  There  is  a  pronounced  rebound  in  FP 
amplitude  when  the  stimulus  intensity  is  returned  to  the  probe  intensity.  E)  A  simple  dynamic 
synapse  model  incorporating  recruitment  of  presynaptic  inputs  is  able  to  reproduce  the  effects  of 
DCS  on  synaptic  efficacy  and  the  post-DCS  rebound  effect. 
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Figure  11.  DCS  modulates  synaptic  efficacy  post-adaptation  through  somatic  polarization  alone.  A) 
Schematic  of  orthodromic  stimulation  of  L23  during  DCS.  Current  flow  along  the  somatodendritic 
axis  and  perpendicular  to  the  horizontally  oriented  axons  selectively  polarizes  the  postsynaptic 
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cells  and  not  the  axons.  B)  Current  flow  along  the  somatodendritic  axis  during  orthodromic 
stimulation  of  the  horizontal  L2/3  pathway  modulates  synaptic  efficacy  Independent  of  axonal 
recruitment.  C)  Anodal  DCS  facilitates  synaptic  efficacy  and  cathodal  DCS  depresses  efficacy. 


A  change  In  the  fiber  volley  (FV)  amplitude,  which  Is  a  compound  action  potential  that  reflects  the  number 
of  activated  axons,  provides  direct  evidence  for  recruitment  with  DCS.  EFs  along  the  ascending  fibers  In 
motor  cortex  (L5^2/3)  decreased  FV  amplitude  by  -15.5±10%  with  -20  V/m  (n=30)  and  Increased  FV 
amplitude  by  13.9±7.1%  with  +20  V/m  (n=30)  (Figure  20A).  Similarly,  In  hippocampal  CA1  area,  EFs 
along  Schaffer  collaterals  Increased  FV  amplitude  by  10.6±10.4%  (p<0.05)  with  -20  V/m  (n=20)  and 
decreased  FV  amplitude  by  -5.5±8%  (p<0.05)  with  +20  V/m  (n=20)  (Figure  20B).  The  polarity  of 
modulation  during  DCS  reflects  axonal  polarization  such  that  depolarized  regions  of  the  axon  (proximal  to 
the  negative  electrode)  correspond  to  an  Increase  In  FV  amplitude.  In  motor  cortex,  there  Is  a  direct 
positive  correlation  (r:  0.92,  Cl:  [0.88,  0.95])  between  change  In  FP  amplitude  and  change  In  FV 
amplitude.  Recruitment  of  afferent  axons  during  DCS  In  brain  slices  may  contribute  to  the  modulation  of 
synaptic  activity  (Figure  20C). 


A 


+ 


la 

o 

o 

o 

LU 


Motor  cortex 


C  Motor  cortex 

±  20  V/m 


45 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Figure  12.  Current  flow  along  axons  during  DCS  modulates  fiber  volley  in  rat  primary  motor  cortex 
and  hippocampus.  A)  Current  flow  along  the  somatodendritic  axis  during  orthodromic  stimulation 
of  the  ascending  L5-*2/3  pathway  polarizes  afferent  axons  and  modulates  axonal  recruitment 
during  DCS.  The  amplitude  of  the  fiber  volley  (FV)  reflects  the  number  of  axons  firing  action 
potentials  during  orthodromic  stimulation.  B)  Similarly,  current  flow  along  the  afferent  axons  of  the 
Schaffer  collateral  pathway  modulates  the  fiber  volley  amplitude.  C)  The  change  in  fiber  volley 
amplitudes  is  directly  correlated  with  the  change  in  postsynaptic  field  potential  amplitude  in  the  rat 
primary  motor  cortex.  These  results  show  polarizing  afferent  axons  changes  the  number  of  active 
presynaptic  inputs  and  directly  modulates  synaptic  current. 


Discussion 

Our  results  show  that  DCS  produces  a  sustained  modulation  of  synaptic  efficacy  that  is  cumulative  over 
the  duration  of  stimulation  when  applied  during  ongoing  synaptic  activity.  Ongoing  synaptic  activity 
produces  well  established  changes  in  synaptic  dynamics  (Abbott  et  al.,  1997;  Varela  et  al.,  1997;  Zucker 
and  Regehr,  2002;  Abbott  and  Regehr,  2004;  Richardson  et  al.,  2005),  with  depression  dominating  in 
motor  cortex.  We  addressed  the  question  if/how  synaptic  dynamics  would  affect  modulation  by  DCS 
(Figure  14).  We  show  synaptic  depression  contributes  to  moderate  loss  in  acute  sensitivity  to  DCS  (from 
1.1%  per  V/m  to  0.55%  per  V/m;  Figure  15E).  Flowever,  there  is  a  significant  modulation  that  is  robust 
across  stimulation  frequency  (Figure  16)  and  extends  to  naturalistic  input  (Figure  17)  and  post-adaptation 
states  (Figure  18).  Importantly,  a  cumulative  gain  in  synaptic  potentials  amplifies  the  DCS  effect  because 
the  modulation  of  synaptic  efficacy  by  DCS  is  sustained  (Figure  15F).  Ongoing  synaptic  activity  combined 
with  sustained  DCS  is  an  important  feature  of  electrical  stimulation  -  quantified  with  a  model  of  synaptic 
dynamics  (Figure  17,  Figure  18e).  We  also  demonstrate  that  afferent  axon  polarization  can  drive  the 
changes  in  synaptic  activity  during  DCS  (Figure  20),  adjunct  to  traditional  somatic  polarization  (Figure  19) 
and  axon  terminal  mediated  mechanisms  (Flause,  1975;  Rahman  et  al.,  2013).  A  targeted  brain  region 
may  be  influenced  by  both  direct  modulation  of  local  neuronal  processes  and  by  polarizing  afferent  brain 
regions.  Although  our  experiments  are  in  vitro,  these  results  have  important  consequences  for  how  the 
outcomes  of  diffuse  and  weak  electrical  stimulation  are  interpreted,  and  suggest  these  features  can  be 
advantageous. 


46 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


tDCS  using  conventional  montages  with  two  large  pad  electrodes  produces  diffuse  current  flow  through 
the  cortex  and  subcortical  brain  regions  (Datta  et  al.,  2009;  Bikson  et  al.,  2010).  For  example,  the 
common  M1-SO  tDCS  montage  (Datta  et  al.,  2009)  produces  diffuse  current  flow  not  only  on  the  nominal 
brain  target,  but  also  in  premotor  and  supplementary  motor  cortex  (Figure  21A),  in  addition  to  other 
cortical  and  subcortical  brain  regions  (Boizoni  et  al.,  2013a;  Boizoni  et  al.,  2013b;  Baczyk  and  Jankowska, 
2014).  Consequently,  the  cellular  effects  of  tDCS  are  not  localized  to  just  the  targeted  brain  region  but 
can  extend  to  upstream  brain  regions.  In  this  way  diffuse  tDCS  may  recruit  a  distributed  brain  network,  for 
instance  the  network  of  brain  regions  involved  in  motor  planning  or  execution.  Flowever,  diffusivity  in  itself 
is  not  a  substrate  for  producing  specific  behavioral  or  clinical  changes,  and  needs  to  be  considered  in  the 
context  of  endogenous  activity  between  connected  brain  regions.  Modulation  of  upstream  brain  regions 
can  influence  the  cooperative  action  of  neurons  and  thus  directly  affect  downstream  processes. 

The  notion  of  diffuse  and  sustained  current  flow  being  advantageous  features  of  tDCS  further  extends  to 
how  the  weak  electric  fields  produced  by  tDCS  are  amplified  by  brain  activity.  The  difference  between 
resting  membrane  potential  and  spike  threshold  for  a  typical  cortical  pyramidal  cell  can  be  15-20  mV. 
DCS,  however,  results  in  a  small  change  in  membrane  potential  (-0.2  mV  per  V/m)  at  the  soma. 
Amplifying  synaptic  activity  through  cooperative  inputs  during  DCS  provides  a  plausible  mechanism  by 
which  weak  EFs  can  significantly  modulate  the  postsynaptic  response  (Equation  1).  An  increase  in 
presynaptic  firing  rate  or  the  synchronous  activation  of  presynaptic  cells  during  diffuse  current  flow  can 
directly  increase  postsynaptic  depolarization  and  firing  probability  (Figure  21  A,  21 B). 

To  formalize  our  hypothesis  that  diffuse  current  flow  may  be  advantageous  for  neuronal  information 
processing,  we  present  a  statistical  theory  (Equation  5)  for  increasing  synaptic  strength  through 
coincident  pre-  and  postsynaptic  spikes.  We  consider  the  number  of  presynaptic  spikes  (iVp^e)  the 
number  of  postsynaptic  spikes  (Wpo^t)  in  discrete  time  bins  (iV^ins).  The  hypergeometric  distribution 
estimates  the  coincidence  count  probability  P(k\Npre,Npost>Nbins)  of  observing  k  coincident  pre-  and 
postsynaptic  spikes  from  a  given  number  of  bins  and  having  observed  Mpg^t  postsynaptic  spikes  and  Np^e 
presynaptic  spikes.  By  varying  the  number  of  presynaptic  spikes  or  postsynaptic  spikes  we  have  an 
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estimate  of  the  number  of  coincidences.  The  probability  distribution  of  coincidences  closely  matches  the 
true  number  of  coincidences  in  simulations  of  neuronal  activity  (Figure  21 C,  dotted  vertical  line  is 
estimated  number  of  coincidences  from  simulations  of  neuronal  activity).  During  tDCS  a  change  in 
presynaptic  firing  rate  (as  a  result  of  diffuse  current  flow),  or  postsynaptic  firing  rate,  or  both,  can  increase 
the  coincidence  probability  and  thus  promote  synaptic  strengthening  between  connected  and  active  brain 
regions.  While  weak  electric  fields  cannot  induce  firing  on  their  own  they  exert  a  modulatory  effect  on 
active  networks  and  increase  the  probability  of  synaptic  strengthening. 
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Figure  13.  tDCS  produces  diffuse  current  flow  in  the  brain  which  increases  the  probability  of 
coincident  pre-  and  postsynaptic  inputs.  A)  The  electric  field  magnitude  is  plotted  in  a  finite 
element  model  of  tDCS  with  anode  over  the  primary  motor  cortex  (Ml)  and  cathode  over 
supraorbital  region  (SO).  The  electric  field  magnitude  distribution  in  brain  regions  adjacent  and 
upstream  to  Ml  is  comparable  to  or  higher  than  in  Ml .  B)  The  primary  motor  cortex  is  synaptically 
connected  to  the  supplementary  motor  area  (SMA)  and  premotor  cortex  (PMC),  which  drives 
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motor  activity.  The  thickness  of  the  arrows  indicates  the  relative  connection  strength  from  fMRI 
studies  of  this  functional  network.  Modulating  presynaptic  activity  in  these  upstream  brain  regions 
during  tDCS  may  drive  postsynaptic  activity  in  M1.  C)  DCS-induced  membrane  polarization 
changes  the  likelihood  of  firing.  An  increase  in  the  presynaptic  activity  (or  postsynaptic  activity) 
increases  the  probability  of  coincident  pre-  and  postsynaptic  action  potentials.  The  vertical  grey 
dotted-line  indicates  the  simulated  number  of  coincidences  in  a  model  l&F  neuron  receiving 
Poisson  distributed  synaptic  inputs.  The  hypergeometric  probability  distribution  can  be  used  to 
approximate  the  estimated  number  of  coincidences.  Combined,  these  results  suggest  the 
probability  of  coincident  inputs  can  be  directly  modulated  when  presynaptic  or  postsynaptic 
activity  is  modulated  by  diffuse  current  flow. 


Implications  and  limitations  of  the  theory  of  diffuse  actions  of  tES 

Future  studies  in  both  humans  and  animals  will  have  to  address  the  role  of  upstream  brain  regions  in 
regulating  the  function  and  output  of  the  target  brain  region  to  test  our  hypothesis.  Notable,  using 
convention  tDCS  montages  in  humans,  peak  brain  current  densities  are  often  between,  rather  than  under 
electrodes;  meaning  significant  diffuse  neuromodulation  is  expected.  Our  statistical  theory  predicts 
changes  in  presynaptic  firing  rate  will  increase  the  coincidence  of  pre-  and  post-synaptic  APs  and 
therefore  increase  the  likelihood  of  plasticity.  Finally,  our  experimental  results  can  provide  some  insights 
into  the  key  features  of  a  typical  tDCS  protocol.  tDCS  produces  diffuse  current  flow  in  the  brain,  needs  to 
be  on  for  an  extended  period  of  time  to  integrate  synaptic  activity,  and  may  have  to  be  matched  with  a 
task.  Our  results  suggest  these  features  may  promote  synaptic  plasticity  because,  as  we  demonstrate, 
changes  in  synaptic  efficacy  are  sustained  for  the  duration  of  DCS  and  helps  synaptic  inputs  integrate. 
Additionally,  diffuse  current  flow  may  enable  recruitment  of  afferent  inputs  in  an  active  functional  network. 
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Section  3:  Electrical  stimulation  accelerates  and  boosts  the  capacity  for  synaptic 

learning 


Introduction 

The  functional  specificity  of  local  synaptic  connections  is  fundamental  to  information  processing  and 
learning  in  the  brain.  Increased  connectivity  between  neurons  with  similar  functional  selectivity  has  been 
reported  across  brain  regions.  Moreover,  synaptic  strengths  within  and  between  these  selective 
populations  are  experience  dependent  and  modifiable  during  a  training  task.  Noninvasive  electrical  brain 
stimulation  may  have  task-specific  effects  where  learning  in  a  functionally  selective  population  benefits 
from  a  combination  of  training  and  the  increased  neuronal  excitability  during  stimulation.  How 
endogenous  synaptic  plasticity  in  local  circuits  is  modulated  by  electrical  stimulation  is  not  well 
understood,  especially  in  the  context  of  task-specific  learning. 

We  found  that  electric  fields  produce  a  voltage-  and  NMDA-dependent  modulation  of  synaptic  strength 
when  applied  during  theta-burst  stimulation  in  rat  hippocampus,  accelerates  synaptic  plasticity  with 
repeated  stimulation,  and  increases  the  ceiling  for  synaptic  learning.  This  mechanism  provides  a  means 
by  which  repeated  stimulation  during  training  can  increase  the  capacity  for  storing  memories  and  gives 
rise  to  the  functional  specificity  of  electrical  brain  stimulation. 


Methods 

All  animal  experiments  were  carried  out  in  accordance  with  guidelines  and  protocols  approved  by  the 
Institutional  Animal  Care  and  Use  Committee  (lACUC)  at  The  City  College  of  New  York,  CUNY  (Protocol 
No:  846.3). 

Brain  slices  including  a  part  of  the  hippocampus  were  prepared  from  male  young  adult  Wistar  rats  aged  3 
to  5  weeks  old,  which  were  deeply  anesthetized  with  ketamine  (7.4  mg/kg)  and  xylazine  (0.7  mg/kg) 
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applied  intraperitoneally  (i.p.)  and  sacrificed  by  cervical  dislocation.  The  brain  was  removed  and 
immersed  in  chilled  (2-6  °C)  artificial  cerebrospinal  fluid  (ACSF)  containing  (in  mM):  126  NaCI,  4.4  KCI, 
1.25  NaH2P04,  1  MgS04,  2  CaCI,  24  NaHCOS  and  10  D-glucose,  bubbled  with  a  mixture  of  95%  02-5% 
C02).  Hippocampal  slices  (400  pm  thick)  were  cut  using  a  vibrating  microtome  and  transferred  to  a 
holding  chamber  for  at  least  30  min  in  ambient  temperature.  Slices  were  then  transferred  to  a  fluid-gas 
interface  chamber  perfused  with  warmed  ACSF  (30.0±0.5  °C)  at  1.9  mL/min.  The  humidified  atmosphere 
over  the  slices  was  saturated  with  a  mixture  of  95%  02-5%  002.  Recordings  started  3-4  h  after 
dissection. 

Orthodromic  stimulation  of  Schaffer  collateral  pathway  within  the  CA1  produced  field  potentials  in  the  LM 
layer.  A  recording  electrode  was  placed  in  the  LM  layer  and  a  bipolar  platinum/stainless  steel  stimulating 
electrode  was  positioned  -500  pm  lateral  to  the  recording  electrode.  Uniform  extracellular  EFs  (±10  and 
20  V/m)  were  generated  by  passing  constant  current  (D/A  driven  analog  follower;  A-M  Systems,  WA, 
USA)  between  two  large  Ag-AgCI  wires  positioned  in  the  bath  across  the  slice  starting  0.5  s  before  theta 
burst  stimulation.  For  all  experiments,  we  used  a  two-tailed  t-test  and  a  bootstrap  hypothesis  test  to  test 
for  significance. 

LTP  was  induced  along  the  Schaffer  collateral  pathway  using  theta  burst  stimulation  (4  pulses  at  100  Hz 
repeated  for  15  bursts  at  5  Hz).  DCS  was  applied  for  3  seconds  during  the  entire  TBS  session.  For 
repeated  stimulation  experiments,  DCS  and  TBS  was  applied  every  40  min  -  1  hr  until  LTP  was 
saturated.  Most  slices  tested  were  saturated  after  3-4  sessions  of  stimulation.  I  also  tested  the  effects  of 
AC  stimulation  when  TBS  is  applied  at  different  phases  of  stimulation.  Each  of  the  15  bursts  in  a  single 
TBS  session  arrived  either  at  the  peak,  trough,  or  at  another  phase  of  AC  stimulation.  In  functional 
specificity  experiments  I  simultaneously  stimulated  both  the  Schaffer  collateral  and  performant  pathways. 
Schaffer  collaterals  were  stimulated  with  TBS  while  the  perforant  pathway  was  stimulated  with  low 
frequency  stimulation  (1 5  pulses  at  5  Hz).  Each  of  the  1 5  pules  arrived  in  the  middle  of  a  burst. 
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Results 


DCS  modulates  endogenous  synaptic  plasticity 

To  mimic  the  learning  during  a  training  task,  we  applied  theta-burst  stimulation  (TBS,  4  pulses  at  100  Hz 
repeated  for  15  bursts  at  5  Hz)  in  the  hippocampal  Schaffer  collateral  pathway  to  induce  synaptic  learning 
through  long-term  potentiation  (LTP).  To  test  the  modulatory  effect  of  electrical  stimulation  we  subjected 
CA1  pyramidal  neurons  to  depolarizing  and  hyperpolarizing  electric  fields  (3  sec;  1,  10,  20  V/m;  Fig.  la). 
We  found  that,  depending  on  the  electric  field  polarity,  the  synaptic  strength  could  be  either  increased  or 
decreased  if  the  field  was  applied  during  TBS.  Soma  depolarizing  direct  current  stimulation  (DCS)  during 
LTP  induction  increased  the  field  EPSP  amplitude  relative  to  baseline  (control:  1.25±0.04,  +20  V/m: 
1.46±0.02),  whereas  soma  hyperpolarizing  DCS  decreased  the  magnitude  of  LTP  (-20  V/m:  1.22±0.01; 
Fig.  1b).  Weak  membrane  polarization  with  DCS  alone  is  not  sufficient  to  induce  plasticity  along  the 
Schaffer  collateral  pathway  in  the  absence  of  synaptic  stimulation. 


DCS  accelerates  and  boosts  the  capacity  for  synaptic  learning 

The  ability  to  modulate  endogenous  synaptic  plasticity  with  electrical  stimulation  provides  a  substrate  for 
enhancing  task-specific  learning  and  prompted  us  to  investigate  whether  repeated  stimulation  sessions 
can  accelerate  synaptic  learning.  Indeed,  we  found  that  repeated  theta  burst  stimulation  (TBS  repeated 
every  60  minutes)  combined  with  DCS  increases  the  rate  of  synaptic  strengthening  until  local  synaptic 
connections  are  saturated  (Fig.  22e).  Inhibiting  synaptic  current  with  the  selective  NMDA  receptor 
antagonist  2-amino-5-phosphonopentanoic  acid  (AP5)  reduced  the  capacity  for  synaptic  learning  and 
abolished  the  modulatory  effect  of  DCS.  Remarkably,  DCS  also  boosted  the  ceiling  of  LTP  (Fig.  22e), 
which  could  be  attributed  to  the  activity-dependent  release  of  extracellular  BDNF  that  has  been  observed 
during  electrical  stimulation  (Fritsch  et  al.,  2010;  Ranieri  et  al.,  2012).  TBS  in  the  presence  of  extracellular 
BDNF  has  also  been  shown  to  increase  the  ceiling  of  synaptic  plasticity. 
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Timing  of  synaptic  inputs  during  ACS  controls  LTP  strength 

The  magnitude  of  LTP  depends  both  on  the  postsynaptic  membrane  potential  and  the  timing  of  synaptic 
inputs  relative  to  oscillations  in  membrane  potential.  A  simple  voltage-based  learning  model  showed  that 
membrane  polarization  controls  the  postsynaptic  firing  probability  leading  to  synaptic  strengthening 
through  coincident  pre-  and  postsynaptic-spikes  (Fig.  22D).  Electrical  stimulation  may  directly  control  the 
magnitude  of  LTP  through  varying  degrees  of  polarization.  In  order  to  study  the  voltage-dependence  of 
LTP  during  electrical  stimulation  we  applied  TBS  at  different  phases  of  theta-frequency  membrane 
potential  oscillation  induced  by  alternating  current  stimulation  (ACS).  We  found  that  the  timing  of  synaptic 
inputs  during  membrane  potential  oscillation  controls  the  magnitude  of  LTP.  TBS  timed  to  the  peak 
depolarization  phase  enhanced  LTP  with  biphasicACS  (figure  22E;  control:  1.25±0.6;  +20  V/m:  1.59±0.03 
relative  fEPSP  slope)  and  with  monophasic  ACS.  Synaptic  inputs  timed  to  the  trough  of  the  membrane 
potential  oscillation,  when  the  cell  is  maximally  hyperpolarized,  decreased  LTP  (Figure  22E;  control: 
1.25±0.6;  -20  V/m:  1.22±0.05  relative  fEPSP  slope).  The  dose-response  (change  in  synaptic  strength  vs 
electric  field  intensity)  has  a  monotonic  and  asymmetric  relationship  (Figure  22F).  These  results  extend 
previous  findings  that  timing  of  synaptic  inputs  during  endogenous  membrane  potential  oscillations 
control  the  magnitude  of  plasticity. 
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Figure  14.  Modulation  of  LTP  with  direct  and  alternating  current  stimulation  accelerates  and  boosts 
the  ceiling  for  synaptic  learning.  A)  Schematic  of  synaptic  inputs  along  the  Schaffer  collateral 
pathway  during  electrical  stimulation  suing  parallel  Ag-AgCI  electrodes.  B)  The  polarity  of  direct 
current  stimulation  (DCS)  controls  the  magnitude  of  plasticity  induced  by  theta-burst  stimulation 
(TBS).  C)  Repeated  TBS  applied  during  DCS  accelerates  synaptic  learning  towards  saturation  of 
synaptic  weights.  Repeated  sessions  of  TBS  with  soma  depolarizing  DCS  (+20  V/m)  increases 
the  ceiling  for  synaptic  learning.  LTP  induction  and  modulation  is  NMDA-dependent.  D)  A  voltage- 
based  learning  rule  in  a  spiking  neuron  model  receiving  correlated  Poisson-distributed  synaptic 
inputs  to  induce  synaptic  learning  is  able  to  reproduce  the  asymmetric  and  voltage-dependent 
potentiation  observed  in-vitro.  DCS  increases  the  spiking  probability  and  promotes  plasticity.  E) 
TBS  is  applied  at  the  peak  or  trough  of  membrane  potential  oscillations  induced  by  alternating 
current  stimulation  (ACS).  Synaptic  inputs  at  the  peak  of  membrane  polarization  boosts  LTP  and 
inputs  at  the  trough  decrease  LTP.  F)  Timing  of  TBS  applied  at  different  phases  of  membrane 
potential  oscillation  controls  the  magnitude  of  LTP. 


Functional  specificity  and  activity-dependent  modulation 

Transcranial  electrical  stimulation  induces  a  diffuse  electric  field  in  the  brain  and  simultaneously  affects 
multiple  brain  regions  during  a  training  task.  TES  induces  significant  and  similar  magnitude  electric  fields 
across  multiple  brain  regions  using  the  traditional  electrode  montage  to  target  primary  motor  cortex.  Flow 
can  a  diffuse  electric  field  have  functionally  specific  effects?  We  hypothesized  the  functional  specificity  of 
TES  arises  from  how  electric  fields  affect  endogenous  synaptic  plasticity.  Active  synaptic  pathways  during 
a  training  task  may  be  modulated  more  than  inactive  or  weakly  active  pathways. 


We  adapted  our  model  to  test  the  functional  specificity  of  electrical  stimulation  by  simulating  a  single 
postsynaptic  neuron  receiving  simultaneous  input  signals  from  two  independent  synaptic  pathways.  The 
strong  pathway  was  activated  with  TBS  to  produce  strong  synaptic  plasticity.  The  weak  pathway  was 
activated  at  a  low-frequency  (5  FIz)  applied  in-phase  with  the  TBS.  Synaptic  inputs  along  the  weak 
pathway  arrive  during  a  critical  period  for  LTP  induction  where  the  TBS  has  already  depolarized  the 
postsynaptic  cell  and  inputs  at  a  low-frequency,  which  typically  induce  depression,  are  now  able  to  evoke 
a  postsynaptic  spike  and  LTP.  Linder  control  conditions  both  pathways  showed  increased  synaptic 
strength  after  training  with  simultaneous  TBS  and  LFS  >  zlw^®®troi)-  Electrical  stimulation  during 

training  modulated  plasticity  in  both  pathways  but  had  a  larger  modulatory  effect  on  the  pathway  with 

stronger  synaptic  plasticity  ( — ^ — ).  Active  pathways  undergoing  strong  synaptic  plasticity 

^'^control  ^'^control 
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may  be  preferentially  modulated  by  electrical  stimulation  during  a  training  task  -  leading  to  a 
discrimination  between  active  and  inactive  functional  networks  during  diffuse  electrical  stimulation. 

We  evaluated  the  modeling  results  experimentally  by  reproducing  the  stimulation  protocol  in-vitro.  A 
conditioning  stimulus  was  applied  to  the  Schaffer  collateral  fibers;  the  test  stimulus  was  applied  to  a 
subicular  input  on  the  opposite  side  of  the  recording  site  (Figure  23A).  Each  shock  of  the  test  input  was 
superimposed  on  the  middle  of  each  burst  of  the  conditioning  input.  The  low  frequency  test  stimulus  was 
first  applied  alone  and  did  not  produce  long-lasting  changes  by  itself  or  in  the  presence  of  an  electric  field 
(Figure  23B).  The  conditioning  site  was  then  stimulated  alone,  which  elicited  homosynaptic  LTP  of  the 
conditioning  pathway  but  did  not  alter  amplitude  of  responses  to  the  test  input  (Figure  23C).  Stimulating 
the  conditioning  site  alone  during  DCS  modulated  synaptic  strength  in  the  conditioning  pathway  but  not  in 
the  inactive  test  pathway  (Figure  24D). 

When  test  and  conditioning  inputs  were  activated  in  phase,  the  conditioning  pathway  was  potentiated  as 
before  and  there  was  an  additional  associative  potentiation  of  the  weakly  activated  test  input  synapses 
(conditioning:  1.21±0.06,  test:  1.16±0.09).  Soma  depolarizing  DCS  boosted  synaptic  strength  in  both 
conditioning  and  test  pathways  while  the  two  pathways  were  simultaneously  activated  (conditioning: 
1.81±0.25,  test:  1.45±0.02).  Electrical  stimulation  had  a  greater  effect  in  modulating  synaptic  plasticity  in 
the  conditioning  pathway  than  the  weakly  activated  test  pathway  -  as  predicted  by  the  model.  The 
preferential  modulation  of  active  local  inputs  is  especially  evident  when  comparing  the  change  in  synaptic 
weights  during  DCS  against  control  Aw^g^troi)  foi"  s^ch  pathway. 
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Figure  15.  Electrical  stimulation  differentially  modulates  strong  and  weak  synaptic  plasticity.  A) 
Model  neuron  receiving  simultaneous  inputs  from  two  independent  signal  groups.  TBS  along  the 
strong  pathway  induces  strong  synaptic  plasticity  and  conditions  a  weak  pathway  to  potentiate 
during  LFS.  B)  Model  prediction  of  the  differential  modulation  of  strong  and  weak  synaptic 
plasticity  during  simultaneous  synaptic  plasticity.  C)  Schematic  of  experimental  setup  where  TBS 
is  applied  to  Schaffer  collaterals  and  LFS  to  the  temporo-ammonic  pathway.  D)  Summary  of 
experimental  results.  Electrical  stimulation  has  a  larger  modulatory  effect  on  LTP  in  strong 
pathways  (conditioning  input)  than  on  weak  pathways  (test  input).  E)  Summary  of  changes  in 
synaptic  strength  for  each  signal  group  during  DCS  and  in  control  conditions  (no  electrical 
stimulation).  F)  Changes  in  synaptic  strength  for  each  signal  group  under  all  tested  conditions. 


Flow  weak  currents  accelerate  task  specific  learning  in  humans  is  still  not  well  understood.  We  found  that 
electric  fields  produce  a  voltage-  and  NMDA-dependent  modulation  of  synaptic  strength  when  applied 
during  theta-burst  stimulation  in  rat  hippocampus,  accelerates  synaptic  plasticity  with  repeated 
stimulation,  and  increases  the  ceiling  for  synaptic  learning.  This  mechanism  provides  a  means  by  which 
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repeated  stimulation  during  training  can  increase  the  capacity  for  storing  memories  and  gives  rise  to  the 
functional  specificity  of  electrical  brain  stimulation. 
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Section  4:  Modulation  of  synaptic  plasticity  with  direct  current  stimulation  is 
activity-dependent  and  differs  with  dendritic  location. 


Introduction 

Transcranial  direct  current  stimulation  (tDCS)  applies  a  weak  constant  current  of  2mA  or  less  across  the 
scalp.  This  apparently  simple  technique  is  currently  under  investigation  for  a  wide  variety  of  conditions, 
including  psychiatric  disorders,  neurorehabilitation  and  cognitive  enhancement  (Coffman  et  al.,  2014; 
Gomez  Palacio  Schjetnan  et  al.,  2013;  Kuo  et  al.,  2014).  Stimulation  is  often  paired  with  a  training  task, 
leading  to  task-specific  enhancements  in  learning  performance  (Coffman  et  al.,  2014;  Gill  et  al.,  2015). 
Despite  the  observation  of  pharmacological,  neuro-physiological  and  imaging  effects  in  humans  (Brunoni 
et  al.,  2012)  and  animals  (Pelletier  and  Cicchetti,  2015),  a  coherent  picture  of  the  relevant  cellular 
mechanisms  is  yet  to  emerge. 

Long-term  learning  and  memory  are  thought  to  be  mediated  by  synaptic  plasticity  (Takeuchi  et  al.,  2014) 
and  training  paradigms  in  humans  presumably  influence  learning  by  inducing  synaptic  plasticity 
(Draganski  and  May,  2008).  Despite  the  common  practice  of  applying  tDCS  during  training,  cellular 
effects  of  DCS  applied  during  endogenous  synaptic  plasticity  induction  remain  largely  unexplored. 
Instead,  the  majority  of  research  has  analyzed  effects  when  DCS  precedes  plasticity  induction  (Podda  et 
al.,  2016;  Ranieri  et  al.,  2012;  Rohan  et  al.,  2015),  or  is  paired  with  endogenous  activity  otherwise  not 
known  to  induce  plasticity  (Fritsch  et  al.,  2010;  Marquez-Ruiz  et  al.,  2012;  Monai  et  al.,  2016).  Here  we 
are  interested  in  the  effects  of  DCS  applied  during  training,  i.e.  concurrent  with  synaptic  plasticity 
induction.  As  a  model  of  endogenous  synaptic  plasticity,  we  induced  long-term  potentiation  (LTP)  and 
depression  (LTD)  using  canonical  induction  protocols  (pulse  trains  delivered  to  Schaffer  collateral 
synapses  in  CA1  of  rat  hippocampal  slices).  By  sweeping  across  induction  frequencies  we  capture  a 
frequency-response  function  (FRF),  which  has  been  widely  used  to  study  the  LTP/LTD  threshold  and  the 
predictions  of  the  Bienenstock,  Cooper  and  Munro  (BCM)  theory  of  synaptic  plasticity  (Cooper  and  Bear, 
2012). 
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A  prevailing  mechanistic  explanation  is  that  tDCS  produces  shifts  in  cortical  excitability,  with  anodal 
stimulation  increasing  excitability  and  cathodal  stimulation  decreasing  excitability  (Brunoni  et  al.,  2012). 
This  cortical  excitability  hypothesis  is  rooted  in  physiological  evidence  that  DCS  modulates  membrane 
potential  at  neuronal  somas,  leading  to  changes  in  firing  rate  and  timing  (Bikson  et  al.,  2004;  Purpura  and 
Mcmurtry,  1965;  Radman  et  al.,  2007,  2009;  Reato  et  al.,  2010).  Based  on  these  observations, 
anodal/cathodal  tDCS  are  often  assumed  to  produce  LTP/LTD-like  effects  for  an  entire  brain  region 
(Dayan  et  al.,  2013;  Monte-Silva  et  al.,  2013;  Ridding  and  Ziemann,  2010;  Stagg  and  Nitsche,  2011). 
However,  this  reasoning  ignores  the  complex  gradient  of  membrane  polarization  induced  in  any  neuron 
during  DCS.  For  example,  anodal  stimulation  depolarizes  pyramidal  neuron  somas  but  will  actually 
hyperpolarize  their  apical  dendrites  (Andreasen  and  Nedergaard,  1996;  Bikson  et  al.,  2004).  Given  the 
dependence  of  synaptic  plasticity  on  dendritic  membrane  potential  (Malenka  and  Bear,  2004),  we 
hypothesized  that  DCS  effects  would  track  the  induced  dendritic  polarization  and  vary  with  dendritic 
location.  Consistent  with  this,  we  find  that  cathodal  DCS  facilitates  LTP  in  apical  dendrites  only,  while 
anodal  DCS  facilitates  LTP  in  basal  dendrites  only.  Moreover,  we  find  that  both  anodal  and  cathodal  DCS 
reduce  LTD  in  apical  dendrites.  Importantly,  DCS  did  not  induce  synaptic  plasticity,  but  rather  required  an 
endogenous  source  of  NMDAR  plasticity  to  modulate.  For  the  first  time  to  our  knowledge,  we  directly 
examine  effects  on  synaptic  plasticity  when  DCS  is  applied  during  plasticity  induction.  We  show  that 
these  effects  are  highly  variable  within  a  single  brain  region,  with  outcomes  dependent  on  the  dendritic 
location  and  rate  of  endogenous  synaptic  activity. 

Methods 

All  animal  experiments  were  carried  out  in  accordance  with  guidelines  and  protocols  approved  by  the 
Institutional  Animal  Care  and  Use  Committee  (lACUC)  at  The  City  College  of  New  York,  CUNY  (Protocol 
No:  846.3). 

Hippocampal  brain  slices  were  prepared  from  male  Wistar  rats  aged  3-5  weeks  old,  which  were  deeply 
anaesthetized  with  ketamine  (7.4  mg  kg-1)  and  xylazine  (0.7  mg  kg-1)  applied  I.P,  and  killed  by  cervical 
dislocation.  The  brain  was  quickly  removed  and  immersed  in  chilled  (2-6°C)  artificial  cerebrospinal  fluid 
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(ACSF)  containing  (in  mM):  NaCI,  125;  KCI,  4.4;  NaH2P04,  1;  MgS04,  1.5;  CaCI,  2.5;  NaHCOS,  26;  D- 
glucose,  10;  bubbled  with  a  mixture  of  95%  02-5%  C02.  Transverse  slices  (400  gm  thick)  were  cut 
using  a  vibrating  microtome  (Campden  Instruments)  and  transferred  to  a  holding  chamber  for  at  least  1  h 
at  ambient  temperature.  Slices  were  then  transferred  to  a  fluid-gas  interface  chamber  (Harvard 
Apparatus)  perfused  with  warmed  ACSF  (30.0  ±  0.1  °C)  at  1.0  ml  min-1.  The  humidified  atmosphere  over 
the  slices  was  saturated  with  a  mixture  of  95%  02-5%  002.  Recordings  started  2-3  h  after  dissection. 
Field  excitatory  postsynaptic  potentials  (fEPSPs)  were  evoked  using  a  platinum-iridium  bipolar  stimulating 
electrode  placed  in  either  stratum  radiatum  or  stratum  oriens  of  CA1.  Recording  electrodes  made  from 
glass  micropipettes  pulled  by  a  Sutter  Instruments  P-97  and  filled  with  ACSF  (resistance  1-8  MQ),  were 
placed  in  either  stratum  radiatum  or  stratum  oriens  approximately  250  pm  from  the  stimulating  electrode 
in  CA1  to  record  fEPSPs.  fEPSPs  were  quantified  by  the  average  initial  slope,  taken  during  the  first  0.5 
milliseconds  after  the  onset  of  the  fEPSP.  Stimulus  intensity  was  set  to  evoke  fEPSPs  with  40%  of  the 
maximum  slope,  which  was  determined  at  the  onset  of  recording.  Stable  baseline  fEPSPs  were  recorded 
every  minute  for  at  least  20  minutes  before  any  plasticity  induction  was  applied.  fEPSPs  were  then 
recorded  again  every  minute  for  60  minutes  after  plasticity  induction. 

DCS  was  applied  between  two  parallel  Ag-AgCI  wires  (1  mm  in  diameter  and  12  mm  in  length) 
placed  in  the  bath  on  opposite  sides  of  the  brain  slice  separated  by  10  mm  with  the  recording  site 
approximately  equidistant  from  each  wire.  DCS  wires  were  connected  to  a  current-controlled  analog 
stimulus  isolator  (A-M  Systems)  that  was  controlled  by  PowerLab  hardware  and  LabChart  software  (AD 
Instruments).  Slices  were  oriented  such  that  the  somato-dendritic  axis  of  CA1  pyramidal  neurons  was 
parallel  to  the  electric  field  between  the  DCS  wires  (Figure  1A).  Before  each  recording,  DCS  current 
intensity  was  calibrated  to  produce  a  20  V/m  electric  field  across  each  slice  (typically  100  -  200  pA)  by 
adjusting  the  current  so  that  two  recording  electrodes  separated  by  0.8  mm  in  the  slice  measured  a 
voltage  difference  of  16  mV  (16  mV/0.8  mm  =  20  V/m). 

For  NMDAR  antagonist  experiments,  100  pM  MK-801  (Sigma  Aldrich)  was  included  in  the  ACSF 
perfused  in  the  recording  chamber  throughout  the  experiment.  Because  MK-801  is  an  open  channel 
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blocker,  baseline  fEPSPs  were  recorded  for  at  least  40  minutes  to  ensure  complete  blockade  of  NMDAR 
channels  (Nabavi  et  al.,  201 3a). 

Data  acquisition  and  stimulation  waveforms  were  controlled  with  PowerLab  hardware  and 
LabChart  software  (AD  Systems).  Extracellular  fEPSPs  were  amplified  (lOOx),  low  pass  filtered  (3  kHz), 
and  digitized  (10  kHz).  Synaptic  plasticity  was  quantified  for  each  slice  by  taking  the  average  of  the  last 
ten  fEPSP  slopes  (51-60  minutes  after  plasticity  induction)  and  normalizing  to  the  average  of  baseline 
fEPSP  slopes  (20-1  minutes  before  plasticity  induction).  Two-tailed  student's  t-tests  were  used  to 
determine  the  significance  of  differences  in  synaptic  plasticity  between  conditions,  with  p  <  0.05 
considered  significant. 

Here  we  name  the  polarity  of  stimulation  based  on  the  orientation  of  DCS  relative  to  pyramidal 
neurons.  Following  convention  in  human  tDCS,  DCS  with  the  anode  closer  to  CA1  apical  dendrites  is 
referred  to  as  anodal  stimulation.  Conversely,  DCS  with  the  cathode  closer  to  CA1  apical  dendrites  is 
referred  to  cathodal  stimulation.  Importantly,  apical  dendrites  are  polarized  oppositely  from  basal 
dendrites  and  somas,  regardless  of  DCS  polarity  (Andreasen  and  Nedergaard,  1996;  Bikson  et  al.,  2004; 
Rahman  et  al.,  2013).  So  anodal  DCS  will  depolarize  somas  and  basal  dendrites,  while  hyperpolarizing 
apical  dendrites.  Conversely,  cathodal  DCS  will  hyperpolarize  somas  and  basal  dendrites,  while 
depolarizing  apical  dendrites  (Figure  1A). 

Acute  effects  were  determined  based  on  the  first  response  (two  responses  for  paired  pulse  data) 
during  DCS  and  were  normalized  to  the  average  of  baseline  responses.  Fiber  volley  amplitude  was  taken 
as  the  difference  between  the  trough  of  the  fiber  volley  and  the  mean  of  the  two  surrounding  peaks. 
Paired  pulse  ratio  was  taken  as  ratio  of  the  second  and  first  fEPSP  slopes  during  20  Hz  HFS  (50  ms  inter¬ 
pulse  interval)  in  each  condition. 
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Frequency-response  functions  in  figure  1B  were  each  fit  with  a  double  exponential  function  with 
five  parameters:  S=  c_1  e''(c_2  x)+  c_3  e''(c_4  x)+c_5  where  x  is  the  induction  frequency,  S  is  the 
resulting  synaptic  plasticity,  and  c_1  ...c_5  are  free  parameters. 


Results 

DCS  shifts  the  frequency-response  function 

Trains  of  synaptic  activity  have  conventionally  been  used  to  induce  synaptic  plasticity  in  hippocampal 
slices  (Malenka  and  Bear,  2004).  As  a  model  of  endogenous  synaptic  plasticity,  trains  of  900  pulses  at 
varying  frequencies  (0.5,  1,  5,  20  Hz)  were  applied  to  the  Schaffer  collateral  pathway  synapsing  on  CA1 
apical  dendrites.  Low  frequency  stimulation  (LFS)  generated  LTD  (0.5  Hz:  84.1  ±2.7%,  p<0.001,  n=10;  1 
Hz:  78.9±2.9%,  p<0.0001 ,  n=9),  while  high  frequency  stimulation  (HFS)  generated  LTP  (20  Hz: 
114. 1  ±2.7%,  p<0.001 ,  n=13),  and  an  intermediate  frequency  marked  the  transition  between  LTD  and  LTP 
(5  Hz:  95.9±3.7%,  p=0.30,  n=9).  The  resulting  FRF  (Figure  IB)  maps  the  degree  of  synaptic  activity 
during  induction  to  the  degree  of  resulting  synaptic  plasticity  and  is  consistent  with  existing  literature 
(Cooper  and  Bear,  2012). 

DCS  was  then  applied  during  plasticity  induction  at  each  frequency.  Our  previous  experiments  with  the 
present  preparation  demonstrate  that  cathodal  DCS  depolarizes  CA1  apical  dendrites  (Figure  1A;  Bikson 
et  al.  2004,  figure  10),  and  was  therefore  expected  to  facilitate  LTP  in  this  dendritic  region  (Malenka  and 
Bear,  2004).  DCS  significantly  attenuated  LTD  induced  by  0.5  Hz  (Figure  1C;  97.5±5.9%,  p  =  0.04,  n  =  8) 
and  1  Hz  LFS  (Figure  1 D;  89.2±3.1  %,  p  =  0.03,  n  =  1 0)  and  enhanced  LTP  induced  by  20  Hz  HFS  (Figure 
IF;  128.8±4.7%,  p  =  0.01,  n  =  14).  The  DCS  effect  at  5  Hz  trended  towards  significance  (Figure  IE; 
102.3±2.2%,  p  =  0.14,  n  =  11),  consistent  with  smaller  effects  observed  previously  at  the  threshold 
between  LTP  and  LTD  (Hulme  et  al.,  2012;  Ikegaya  et  al.,  2002;  Lisman,  2001).  The  resulting  DCS  FRF 
was  significantly  shifted  compared  to  control  (F=17.93,  df=1,  p<0.0001).  Fitting  a  double  exponential  to 
the  means  of  each  FRF  highlights  the  DCS-induced  reduction  in  LTP  threshold,  biasing  synapses  toward 
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potentiation  (Figure  1B).  Similar  shifts  of  the  FRF  have  been  associated  with  enhanced  learning  in  cortex 
(Cooper  and  Bear,  2012;  Flulme  et  al.,  2013). 


DCS  effects  on  LTP  depend  on  dendritic  location 

DCS  is  known  to  modulate  the  membrane  potential  of  neuronal  compartments  (Andreasen  and 
Nedergaard,  1996;  Bikson  et  al.,  2004;  Chan  et  al.,  1988;  Radman  et  al.,  2009;  Rahman  et  al.,  2013)  and 
dendritic  membrane  potential  is  known  to  be  a  critical  determinant  of  NMDAR-dependent  plasticity 
(Malenka  and  Bear,  2004).  Other  DCS  effects  in  humans  and  animals  have  been  shown  to  be  NMDAR- 
dependent,  and  it  is  widely  speculated  that  tDCS  exerts  long-term  effects  through  membrane  polarization 
and  NMDARs  (Stagg  and  Nitsche,  2011).  An  important  subtlety  that  is  often  lost  in  this  discussion  is  that 
DCS  will  simultaneously  depolarize  and  hyperpolarize  different  compartments  within  the  same  neuron. 
Indeed,  previous  work  from  our  own  group  with  a  similar  experimental  setup  showed  that  cathodal  DCS 
simultaneously  depolarizes  CA1  apical  dendrites  while  hyperpolarizing  their  basal  dendrites  and  soma. 
Conversely,  anodal  DCS  hyperpolarizes  CA1  apical  dendrites  while  depolarizing  their  basal  dendrites  and 
soma  (Bikson  et  al.,  2004).  We  therefore  expected  that  the  effects  of  anodal  and  cathodal  stimulation 
would  vary  with  dendritic  location.  To  test  this  we  paired  both  anodal  and  cathodal  DCS  with  20  FIz  FIFS 
in  both  CA1  apical  and  basal  dendrites.  In  apical  dendrites,  cathodal  DCS  enhanced  LTP,  while  anodal 
DCS  had  no  significant  effect  (Figure  2C,D;  control:  1 14.1±2.7%,n=13;  cathodal: 
128.8±4.7%,p=0.01  ,n=14;  anodal:  1 1 1 .7±4.5%,p=0.63,  n=8).  In  basal  dendrites,  anodal  DCS  now 
enhanced  LTP  while  cathodal  DCS  had  no  significant  effect  (Figure  2E,F;  control:  148.6±3.6%,n=10; 
cathodal:142.5±5.2%,  p=0.34,n=10;  anodal:  180.4±9.1%,  p<0.01,n=5).  As  expected,  the  effects  of 
anodal  and  cathodal  DCS  were  dependent  on  dendritic  location. 

DCS  effects  are  polarity  dependent  for  LTP  but  not  LTD 

Anodal  and  cathodal  DCS  apply  stimulation  with  opposite  polarity  and  are  canonically  expected  to 
produce  opposite  effects  (Nitsche  et  al.,  2008).  As  reported  above,  we  find  that  cathodal  and  anodal  DCS 
have  asymmetric  effects  on  LTP  for  a  given  dendritic  location.  Moreover,  when  paired  with  IFIz  LFS  we 
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observe  no  polarity  dependence  of  effects.  LTD  is  reduced  by  both  anodal  and  cathodal  DCS,  i.e. 
synaptic  strength  is  increased  compared  to  control  (Figure  2G,H;  control:  78.9±2.9%,n=9;  cathodal: 
89.2±3.1%,p=0.03,n=10;  anodal:  95.6±5.9%,n=8,p=0.04).  These  results  reveal  that  modulation  of 
synaptic  plasticity  by  DCS  depends  on  both  the  physical  location  of  concurrently  active  synapses  (basal 
or  apical  dendrites)  and  the  rate  of  their  activity  (LFS  or  FIFS)  (Figure  2B). 


DCS  effects  require  a  concurrent  endogenous  source  of  NMDAR  plasticity 

tDCS  is  often  applied  under  the  assumption  that  stimulation  alone  is  sufficient  to  induce  plasticity  (Brunoni 
et  al.,  2012;  Pelletier  and  Cicchetti,  2015;  Ridding  and  Ziemann,  2010).  Flowever,  given  the  modest 
effects  on  membrane  potential  for  typical  stimulation  intensities  (Andreasen  and  Nedergaard,  1996; 
Radman  et  al.,  2009;  Rahman  et  al.,  2013),  we  propose  that  DCS  instead  acts  as  a  modulator  of  NMDAR 
plasticity.  DCS  would  therefore  require  a  concurrent  endogenous  source  of  plasticity  to  modulate.  To  test 
this  requirement  we  again  applied  cathodal  DCS,  but  removed  endogenous  NMDAR-dependent  plasticity 
in  two  ways:  first  by  weakening  synaptic  activity  to  well  below  the  plasticity  threshold,  and  second  by 
directly  blocking  NMDAR  current  during  strong  synaptic  activity.  When  applied  during  weak  synaptic 
activity  (30  pulses,  1/60  FIz),  cathodal  DCS  had  no  effect  (Figure  3C;  control:  99.3±1.1%,  n=9;  cathodal 
DCS:  100.8±4.0%,  n=7;  p  =  0.68).  When  paired  with  strong  synaptic  activity  (20  FIz  FIFS)  but  NMDARs 
were  blocked  with  antagonist  MK-801,  cathodal  DCS  also  had  no  effect  (Figure  3B,  control:  92.0±1.6%, 
n=10;  cathodal  DCS:  94.3±2.3%,  n=9;  p  =  0.42).  These  results  suggest  that  DCS  may  act  as  a  modulator 
of  endogenous  synaptic  plasticity,  rather  than  an  inducer  of  de  novo  synaptic  plasticity. 


Acute  effects  of  DCS  on  synaptic  transmission 

To  gain  insight  into  potential  mechanisms  by  which  DCS  may  modulate  synaptic  plasticity,  we  examined 
the  effects  of  DCS  on  several  measures  of  baseline  synaptic  transmission.  Flowever,  one-way  ANOVAs 
yielded  no  significant  effect  of  stimulation  on  fEPSP  slope  (Figure  4A;  F=0.23,df=1  ,p=0.63),  fiber  volley 
amplitude  (Figure  4B,  F=0.33,df=1  ,p=0.57),  or  paired  pulse  ratio  (Figure  4C;  F=0.1 1  ,df=1  ,p=0.74). 
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Discussion 


LTP,  LTD,  and  learning 

There  is  now  strong  evidence  for  a  role  of  both  LTP  and  LTD-like  processes  in  various  types  of  learning 
and  memory  (Broussard  et  al.,  2016;  Connor  and  Wang,  2015;  Dong  et  al.,  2013;  Kuhiman  et  al.,  2014; 
Nabavi  et  al.,  2014;  Rioult-Pedotti,  2000).  At  the  behavioral  level,  learning  is  likely  to  involve  both  of  these 
processes,  with  the  precise  degree  of  each  depending  on  the  specific  behavior.  For  example,  some 
learned  behaviors  directly  require  habituation  to  a  familiar  stimulus  and  are  specifically  dependent  on  LTD 
(Griffiths  et  al.,  2008;  Massey  and  Bashir,  2007).  Other  learned  behaviors  involve  formation  of  new 
associations  and  responses  to  the  environment,  which  require  LTP  and  are  eliminated  by  LTD  (Broussard 
et  al.,  2016;  Nabavi  et  al.,  2014).  We  observed  that  LTP  is  facilitated  in  dendrites  that  are  depolarized  by 
DCS.  This  cellular  DCS  effect  may  contribute  to  enhanced  learning  when  tDCS  is  paired  with  training  that 
induces  plasticity,  such  as  motor  rehabilitation  (Kang  et  al.,  2016;  Rioult-Pedotti,  2000).  Indeed,  similar 
shifts  in  the  FRF  have  been  linked  to  facilitation  of  learning  on  both  theoretical  and  experimental  grounds 
(Cooper  and  Bear,  2012).  We  also  observed  a  reduction  of  LTD  for  both  stimulation  polarities  (Figure 
2B).  One  may  therefore  expect  that  these  effects  would  inhibit  learning  that  requires  LTD. 


Plasticity  dependence  may  underlie  task-specific  effects 

When  tDCS  is  paired  with  training,  the  observed  effects  are  often  specific  to  the  trained  task  (Bikson  et 
al.,  2013;  Gill  et  al.,  2015;  Sipolins  et  al.,  in  review).  While  electrodes  are  typically  placed  over  an 
intended  target  region,  it  is  unlikely  that  task  specificity  is  solely  the  result  of  spatial  selectivity  of  current 
flow.  Even  in  the  most  focal  tDCS  applications  (e.g.  FID-tDCS),  current  flow  through  the  brain  is  diffuse, 
reaching  large  swaths  of  cortex  and  subcortical  structures  (Borckardt  et  al.,  2012;  Datta  et  al.,  2009). 
Moreover,  within  any  particular  brain  region,  there  are  likely  to  be  neurons  involved  in  many  disparate 
memory  engrams  or  behaviors.  The  common  assumption  that  tDCS  induces  plastic  effects 
indiscriminately  (Nitsche  et  al.,  2008),  or  even  at  weakly  active  synapses  (Fritsch  et  al.,  2010),  therefore 
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implies  broad  effects  on  any  cognitive  output  in  the  stimulated  brain  regions.  This  is  at  odds  with  the 
observed  specificity  of  effects.  Instead,  to  explain  task-specificity  tDCS  may  act  as  a  selective  modulator 
of  endogenous  synaptic  plasticity.  Our  results  support  this  hypothesis,  as  DCS  had  no  effect  when 
synaptic  input  was  too  weak  (Figure  3C)  or  when  NMDARs  were  blocked  during  strong  synaptic  input 
(Figure  3B),  indicating  that  synaptic  efficacy  is  modulated  by  DCS  only  when  NMDAR-dependent 
plasticity  is  present.  This  provides  a  basis  for  effects  to  be  task-specific,  as  synapses  associated  with  the 
paired  task  are  more  likely  to  be  undergoing  plasticity  and  therefore  subject  to  modulation  during  tDCS. 
Moreover,  this  predicts  that  tDCS  effects  should  be  enhanced  when  paired  with  tasks  that  induce  synaptic 
plasticity.  Indeed,  there  is  some  evidence  for  this  (Gill  et  al.,  2015;  Martin  et  al.,  2014).  The  precise  role 
of  endogenous  synaptic  activity  in  DCS  effects  remains  an  important  area  for  future  research. 


Dendritic  pciarizaticn  determines  dependence  cn  DCS  pciarity  and 
pathway 

Under  the  conventional  excitability  hypothesis,  ‘anodal  tDCS’  is  assumed  to  produce  inward  cortical 
current  flow,  which  depolarizes  pyramidal  neuron  somas  and  hence  increases  cortical  excitability. 
‘Cathodal  tDCS’  is  soma-hyperpolarizing  and  thus  should  reduce  cortical  excitability  (Nitsche  et  al.,  2008). 
Flowever,  it  is  becoming  increasingly  clear  that  this  reasoning  is  an  oversimplification,  particularly  when  it 
comes  to  long-term  effects  and  learning  (Antal  et  al.,  2007;  Jacobson  et  al.,  2012;  Vines  et  al.,  2008). 
While  effects  on  somatic  membrane  potential  must  still  be  considered,  our  results  here  point  to  a  role  for 
dendritic  membrane  polarization  in  determining  DCS  effects  on  synaptic  plasticity. 

Membrane  polarization  due  to  DCS  can  in  principle  affect  the  function  of  all  voltage-dependent  channels 
distributed  throughout  a  neuron,  particularly  the  relief  of  NMDARs  from  magnesium  blockade.  This 
influence  may  be  most  pronounced  in  dendrites,  where  DCS  has  been  shown  to  modulate  excitability 
involving  multiple  voltage-dependent  channels  (Andreasen  and  Nedergaard,  1996).  While  we  do  not 
directly  measure  membrane  polarization  in  the  present  experiments,  our  group  has  done  this  previously 
with  the  same  preparation,  showing  membrane  polarization  to  be  maximal  in  dendrites  (Bikson  et  al. 
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2004,  figure  10),  with  opposite  polarization  in  apical  and  basal  dendrites.  Indeed  we  observe  modulation 
of  synaptic  plasticity  that  is  consistent  with  this  variable  dendritic,  rather  than  a  singular  somatic 
polarization  effect  (Figure  2). 

Given  that  DCS  effects  can  vary  with  dendritic  location,  tasks  that  activate  pathways  with  different 
dendritic  locations  may  respond  differently  to  the  same  stimulation.  A  lack  of  control  over  the  location  of 
active  pathways  could  therefore  lead  to  highly  variable  results  in  clinical  studies.  Attention  to  dendritic 
polarization  may  therefore  help  to  explain  mixed  effects  observed  in  tDCS  outcomes  (Antal  et  al.,  2007; 
Horvath  et  al.,  2014;  Jacobson  et  al.,  2012;  Saucedo  Marquez  et  al.,  2013).  This  motivates  a  shift  away 
from  the  conventional  focus  on  somatic  polarization,  which  implies  fixed  excitability  changes  for  an  entire 
brain  region. 

While  our  results  are  consistent  with  a  role  for  DCS-induced  dendritic  polarization,  we  cannot  rule  out  the 
involvement  of  other  cellular  DCS  effects,  such  as  on  inhibitory  interneurons,  glia,  neuromodulators 
systems,  or  immune  response  (Pelletier  and  Cicchetti,  2015).  Further  investigation  into  the  involvement  of 
these  systems  is  an  important  area  for  future  work.  The  lack  of  effects  on  fiber  volleys  and  paired  pulse 
ratio  suggest  that  DCS  does  not  affect  recruitment  or  vesicle  release  probability  at  presynaptic  terminals. 
This  is  expected,  as  Schaffer  collateral  fibers  are  oriented  perpendicular  to  the  applied  DCS  electric  field 
vector.  While  previous  studies  have  reported  effects  on  fEPSP  slope,  the  lack  of  an  effect  here  may 
result  from  a  smaller  sample  size  and  weaker  fields  (Bikson  et  al.,  2004;  Kabakov  et  al.,  2012;  Rahman  et 
al.,  2013). 


Low  frequency  stimulation  effects 

The  horizontal  axis  of  the  FRF  is  often  equated  with  the  degree  of  postsynaptic  calcium  influx  during 
induction.  HFS  leads  to  strong  calcium  influx  and  triggers  LTP,  while  LFS  leads  to  moderate  calcium 
influx  and  LTD.  Based  on  this  calcium  control  hypothesis,  we  expected  DCS-induced  dendritic 
polarization  to  modulate  calcium  influx  through  NMDARs  and  produce  horizontal  shifts  in  the  FRF 
(Cooper  and  Bear,  2012).  The  effects  we  observe  with  1  Hz  LFS  may  therefore  be  expected,  as  a 
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horizontal  shift  of  the  FRF  in  either  direction  would  result  in  less  LTD  if  1  FIz  is  near  the  point  of  maximum 
LTD  (minimum  synaptic  strength).  This  interpretation  is  less  adequate  in  accounting  for  the  effect 
observed  at  0.5  FIz  LFS  (Figure  1C),  as  a  left  horizontal  shift  would  produce  more  LTD  at  0.5  FIz. 
Flowever,  recent  evidence  suggests  a  deviation  from  the  calcium  control  hypothesis,  as  LTD  can  be 
induced  by  metabotropic  NMDAR  function  rather  than  calcium  influx  (Aow  et  al.,  2015;  Dore  et  al.,  2015; 
Nabavi  et  al.,  2013a,  2013b;  Stein  et  al.,  2015).  The  calcium  and  voltage  dependence  of  LTD  remains 
controversial  though  (Babiec  et  al.,  2014;  Nabavi  et  al.,  2013b),  making  it  more  difficult  to  interpret  results 
with  LFS.  We  also  note  that  the  duration  of  DCS  was  particularly  long  with  0.5  FIz  LFS  (30  minutes), 
potentially  producing  effects  that  occur  on  longer  time  scales,  such  as  on  protein  synthesis.  For  example, 
priming  of  BDNF  synthesis  at  the  start  of  DCS  (Podda  et  al.,  2016)  may  lead  to  increased  BDNF  release 
later  on  during  DCS,  which  reduces  LTD  (Ikegaya  et  al.,  2002).  Future  experiments  directly  measuring 
calcium  influx  during  these  induction  protocols  may  provide  some  resolution  to  these  issues. 


Effect  asymmetry 

Our  results  demonstrate  an  asymmetric  DCS  effect  on  synaptic  plasticity,  such  that  DCS  was  only  able  to 
increase  synaptic  strength  (enhance  LTP,  reduce  LTD).  Asymmetries  have  been  found  in  other  animal 
studies  (Marquez-Ruiz  et  al.,  2012;  Ranieri  et  al.,  2012)  and  human  studies  (Coffman  et  al.,  2014; 
Jacobson  et  al.,  2012).  In  parallel  work  in  our  lab,  we  find  an  asymmetry  in  acute  DCS  effects  on  cellular 
excitability.  This  nonlinearity  could  be  the  result  of  the  nonlinear  voltage  dependence  of  NMDARs  (Jahr 
and  Stevens,  1990).  Similarly,  these  asymmetries  may  reflect  floor  or  ceiling  effects  of  any  number  of 
cellular  processes,  where  the  endogenous  state  is  such  that  it  can  only  be  modulated  in  one  direction. 


Conclusions  and  context 

DCS  is  likely  to  affect  many  cellular  processes  simultaneously  (Pelletier  and  Cicchetti,  2015).  To  improve 
tDCS  applications  it  will  be  important  to  parse  which  of  these  cellular  effects  are  most  relevant  under 
various  clinical  conditions,  such  that  desirable  outcomes  can  be  reliably  produced.  Previous  studies  in 
animals  (Fritsch  et  al.,  2010;  Marquez-Ruiz  et  al.,  2012;  Monai  et  al.,  2016;  Podda  et  al.,  2016;  Ranieri  et 
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al.,  2012;  Reato  et  al.,  2015;  Rohan  et  al.,  2015)  and  humans  (Liebetanz,  2002;  Nitsche  et  al.,  2003;  Rroji 
et  al.,  2015)  have  implicated  various  effects  related  to  synaptic  plasticity  (NMDAR,  BDNF,  adenosine, 
norepinephrine).  However,  it  remains  unknown  exactly  how  the  DCS  electric  field  vector  interacts  with 
cellular  activity  to  produce  these  effects.  The  brain  slice  preparation  used  here  allows  for  precise  control 
over  the  electric  field  vector  with  respect  to  neuronal  morphology  and  synaptic  activity,  facilitating  a 
bottom-up  approach.  It  will  be  important  for  future  work  to  relate  the  effects  observed  here  to  in  vivo 
tDCS  in  both  rodent  models  and  humans. 

DCS-induced  membrane  polarization  is  often  cited  as  a  source  of  effects  (Stagg  and  Nitsche,  2011),  but 
the  polarization  profile  of  a  neuron  during  DCS  will  be  as  complex  as  its  morphology  (Rahman  et  al., 
2013)  and  the  effects  of  polarization  are  likely  to  vary  with  the  endogenous  activity  of  the  neuron.  We 
demonstrate  here  that  the  effects  of  DCS  vary  with  the  spatial  and  temporal  properties  of  the  endogenous 
synaptic  activity  that  it  is  paired  with  (Figure  2).  These  results  highlight  the  complexity  of  DCS 
interactions  with  neuronal  morphology  and  activity,  motivating  a  departure  from  the  canonical  excitability 
hypothesis,  which  implies  singular  effects  for  an  entire  brain  region. 
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Figure  1.  DCS  shifts  synaptic  plasticity  in  apical  dendrites  towards  potentiation.  A:  Schematic 
depicts  cathodal  DCS  of  a  hippocampal  slice,  with  expected  membrane  polarization  of  CA1  pyramidal 
neuron  (enlarged  at  right;  prediction  based  on  computational  model  as  in  Rahman  et  al.  2013).  Arrow 
indicates  the  direction  of  positive  current  flow  between  electric  field  wires  placed  in  the  recording  chamber 
(horizontal  bars  above  and  below  hippocampal  slice).  B:  Cathodal  DCS  (blue)  shifts  the  frequency- 
response  function  towards  potentiation.  Data  are  fit  with  a  double  exponential  function  to  highlight  the 
change  LTP  threshold.  C-F:  DCS  applied  during  plasticity-inducing  LFS  attenuated  LTD  (C,D)  and 
enhanced  LTP  (F),  but  the  effect  was  not  significant  near  the  crossover  point  between  LTD  and  LTP  (E). 
Sample  fEPSP  traces  are  provided  for  each  condition  (grey/black:  before/after  control;  light  blue/blue: 
before/after  cathodal;  scale  bars:  1  mV,  4  ms).  Synaptic  strength  in  (B)  is  the  average  of  the  last  ten 
normalized  fEPSP  slopes  in  each  condition  (51-60  minutes  post-induction).  Data  for  0.0167  Hz  is  shown 
in  figure  3C.  Grey  bars  indicate  the  duration  of  plasticity  induction  and  concurrent  DCS.  Data  are 
represented  as  mean  ±  SEM  across  slices.  *  =  p  <  0.05. 
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Figure  2.  Mixed  of  effects  of  anodal  and  cathodal  DCS  depend  on  dendritic  location.  A:  Schematic 
depicts  anodal  (red)  and  cathodal  (blue)  DCS  of  a  hippocampal  slice.  Arrows  indicate  the  direction  of 
positive  current  flow  between  DCS  electrodes.  Reconstruction  of  a  CA1  pyramidal  neuron  with  dendritic 
compartments  labeled  at  right.  B:  Modulation  of  synaptic  plasticity  depends  on  polarity,  dendritic  location 
and  rate  of  induction.  C-D:  In  apical  dendrites  cathodal  DCS  enhances  LTP,  but  anodal  has  no 
significant  effect.  E-F:  Changing  dendritic  location  to  basal  dendrites,  anodal  DCS  now  enhances  LTP, 
but  cathodal  DCS  has  no  effect.  G-H:  Changing  induction  frequency  to  1  Hz  LFS  (in  apical  dendrites), 
both  anodal  and  cathodal  reduce  LTD.  C-H:  For  comparison,  DCS  polarity  is  the  same  across  each  row 
of  panels,  while  dendritic  recording  location  is  the  same  within  each  column  of  panels.  Example  traces 
for  each  condition  are  given  in  the  center  of  each  column  (grey/light  blue/pink  traces  are  before  plasticity 
induction;  black/blue/red  traces  are  after  plasticity  induction;  scale  bars:  1  mV,  4  ms).  Insets  depict  the 
reconstructed  CA1  pyramidal  neuron  in  (A)  with  expected  membrane  polarization  (prediction  based  on 
computational  model  as  in  Rahman  et  al.  2013),  induction  and  recording  sites,  and  orientation  of  DCS 
electric  field  vector.  Grey  bars  indicate  the  duration  of  plasticity  induction  and  concurrent  DCS.  Plasticity 
modulation  in  (B)  is  the  resulting  plasticity  in  each  DCS  condition  normalized  to  the  mean  of  the  plasticity 
in  the  corresponding  control  condition.  Data  are  represented  as  mean  ±  SEM  across  slices.  *  =  p  <  0.05. 
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Figure  3.  DCS  requires  an  endogenous  source  of  NMDAR  plasticity.  A:  Schematic  depicts  cathodal 
DCS  (blue)  of  a  hippocampal  slice,  with  expected  membrane  polarization  of  CA1  pyramidal  neuron 
(enlarged  at  right;  prediction  based  on  computational  model  created  as  described  in  Rahman  et  al.  2013). 
B:  Cathodal  DCS  with  20  Hz  HFS  has  no  effect  on  synaptic  strength  when  NMDARs  are  blocked  with 
antagonist  MK-801  C:  Cathodal  DCS  applied  during  synaptic  activity  that  is  too  weak  to  induce  plasticity 
(30  pulses  at  1/60  Hz)  has  no  effect  on  synaptic  strength.  Grey  bars  indicate  duration  of  induction  and 
concurrent  DCS.  Data  are  represented  as  mean  ±  SEM  across  slices. 


Figure  4.  No  significant  effects  on  baseline  synaptic  transmission.  DCS  had  no  significant  effect  on 
fEPSP  slope  (A),  fiber  volley  amplitude  (B),  or  paired  pulse  ratio  (C),  in  apical  or  basal  dendrites,  or  when 
NMDARs  were  blocked  with  MK-801 .  Data  are  represented  as  mean  ±  SEM  across  slices,  ns  =  p  >  0.05. 
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